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SAM Localization of Oscillatory MEG Activities in Human Brain 
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ABSTRACT 

The oscillatory MEG activities have multiple distributed generators in human brain. We have reported the localizations of several oscillatory 
MEG activities by using synthetic aperture magnetometry (SAM). The sources of these oscillatory activities were identified through the permutation 
tests of SAM results from multiple subjects. Source locations of the posterior alpha rhythm were found primarily in the bilateral occipital cortices, 
including the vicinity around the calcarine fissure and in the bilateral parietal cortices, including the vicinity around the parieto-occipital sulcus. We 
suggested that the occipital alpha rhythm has each two clusters at bilateral occipital and parietal cortices. The alpha component of Rolandic mu 
rhythm showed the rebound after stopping finger movements and distributed in the posterior areas of central sulcus whereas the beta and gamma 
components located at the anterior areas of central sulcus. The generators of frontal midline theta rhythm were localized in the medial prefrontal 
cortex and anterior cingulate cortex. We suggest that these areas are related with the focused attentional processing. Current source density 
distribution of sleep spindles was estimated at wide and distributed areas of right frontal and bilateral parietal areas including sensorimotor cortex. 
We speculated that some components of sleep spindles might have common frequency bands and common distribution with Rolandic mu rhythm. 
Knowledge of these oscillatory activities in normal human brain is important for clinical applications to detect the abnormal brain activity, and for 
understanding the psychophysiological meaning of MEG signals. 

KEY WORDS 

Theta rhythm, alpha rhythm, beta rhythm, gamma rhythm, MEG, spatial filter, equivalent current dipole (ECD), SPM, permutation analysis, 
synthetic aperture magnetometry (SAM). 

INTRODUCTION 

Since the first report of Berger in 1929, the alpha rhythm is the well-known oscillatory activity in human EEG, with the frequency around 10Hz, 
which is dominantly distributed at occipital area during eye-closed resting state and attenuate during eye-opening state. The mu and tau rhythms have 
similar frequency bands and distribute dominantly at central and temporal areas respectively. In 1972, Ishihara and Yoshii reported that the theta 
rhythms in EEG were observed at frontal area when the subjects were concentrated on the calculation tasks and they named this activity as Frontal 
Midline Theta Rhythm [Ishihara, 1972], Further studies showed that these phenomenons were observed during the focused attentional states in 
various cognitive tasks. Also in sleep EEG it is possible to observe several oscillatory activities like theta waves in drowsy states, sleep spindles and 
high voltage slow waves in deep sleep. Sleep stages are defined by these EEG phenomenons. 

These oscillatory activities have different frequency bands, scalp distributions and behavior and have been studied in EEG and MEG [Hari, 
1997]. Although it was speculated that these oscillatory activities have wide and spread distributions in human brain, the equivalent current dipole 
(ECD) methods, which localized point-like generators, were widely used for the analysis of the oscillatory activities in MEG studies. We used one of 
the spatial filtering analysis methods, Synthetic Aperture Magnetometry (SAM) [Robinson, 1992] [Ishii, 1999]., and analyzed the steady-state 
oscillatory MEG activity in order to estimate the current source density distributions and visualize these activity on images of anatomical MRI. 
Recently we used SPM99 software (Wellcome Department of Cognitive Neurology, London, UK) for the normalization of SAM images and 
determined the statistical significance of the normalized SAM images by using a non-parametric permutation method [Chau, 2002], 

We analyzed the basic alpha and theta rhythms by using SAM-SPM99-Permutation method. Source locations of the posterior alpha rhythm and 
higher frequency components in beta and gamma bands were found primarily in the bilateral occipital cortices, including the vicinity around the 
calcarine fissure and in the bilateral parietal cortices, including the vicinity around the parieto-occipital sulcus (Figl) [Ishii, 2003]. We suggested that 
the occipital alpha rhythm has two clusters at bilateral occipital and parietal cortices and the beta and gamma components have different distributions 
in occipital areas. The alpha component of Rolandic mu rhythm showed the rebound after stopping finger movements and distributed in the posterior 
areas of the central sulcus, whereas beta and gamma components located at the anterior areas of central sulcus (Figl) [Ishii, 2003]. These results were 
consistent with the results of ECD analysis [Salmelin, 1994], The generators of frontal midline theta rhythm were localized at the medial prefrontal 
cortex and anterior cingulate cortex. Therefore, we suggest that these areas are related with the focused attentional processing [Ishii, 2004], Current 
source density distribution of sleep spindles was estimated at the wide and distributed areas of right frontal and bilateral parietal areas including 
sensorimotor cortex. We speculated that some components of sleep spindles might have common frequency bands and common distribution with 
Rolandic mu rhythm [Ishii, 2003], 
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Pseudo t-valuc 

Figl. SAM-SPM99-Permutation images of the occipital alpha rhythm (left) and the Rolandic mu rhythm (right). 

(n=15, p<0.001) [Ishii, 2003] 


We conclude that analyzing generators of the oscillatory activities are important not only for basic physiological research, but also for clinical 
research to detect the abnormality in the patients with epilepsy, dementia and schizophrenia [Ishii, 2000], Recently, these oscillatory activities were 
investigated in the context of cognitive functions. Since SAM-SPM99-Permutation method can estimate multiple and distributed generators, this 
method has an advantages for these oscillatory activities related to cognitive functions. 

In the past few decades, evoked potential (EP) and event-related potential (ERP) studies were in the main stream of EEG and MEG research and 
oscillatory brain activities have been often treated as “brain noise” and canceled from the recorded signals when applying simple averaging technique 
in the EP and ERP studies. However, the research on oscillatory brain activity is important for clinical applications in order to separate abnormal 
brain activity, and for understanding the psychophysiological meanings of the recorded brain signals. 
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The Neurophysiology of Human Spatial Navigation 

Michael J. Kahana 

University of Pennsylvania, Philadelphia, USA 

The place cells of the rodent hippocampus, and the 3-10 Hz hippocampal theta rhythm, constitute two of the most striking examples of correlations 
between neuronal activity and complex behavior in mammals. Using human intracranial recordings, which can be ethically obtained as part of 
standard neurosurgical evaluations, we examine the neurophysiological basis of human spatial navigation. Field recordings from depth electrodes 
while subjects actively explored a virtual town revealed increased theta activity associated with movement and when searching for objects whose 
locations were not known in both hippocampus and neocortex. Single cell recordings from hippocampal, parahippocampal and frontal regions reveal 
cells in the hippocampus that respond at specific spatial locations and cells in parahippocampal region responding to views of landmarks. Cells also 
responded to the subject’s navigational goals and to conjunctions of place, goal and view. Finally, some cellular responses were phase-locked to the 
locally generated theta rhythm. These data reveal a striking parallel between the neural basis of navigation in human and animals, and demonstrate 
that both map and route type information is used during navigation. 


Multivariate analyses of oscillatory responses of declarative memory and their relationship to event- 

related fields 

E. Diizel 

Department of Neurology II, Otto-von-Guericke University Magdeburg, Germany 

Much of our knowledge about the electrophysiology of declarative memory in humans stems from EEG and MEG signals that were averaged across 
multiple trials resulting in event-related potentials (ERPs) and fields (ERFs). Neural mechanisms of declarative memory in animals, on the other 
hand, have been mostly studied through analyses of oscillatory rhythms. Therefore, the relationship between memory related changes in ERPs/ERFs 
and oscillatory phenomena related to memory is of current interest. It will be shown that the frequency composition of ERF effects of recognition 
memory in humans is in the range of delta and slow theta oscillations. Single trial analyses of phase and amplitude values show that ERF effects of 
recognition memory are related to changes in the phase alignment of delta and slow theta oscillations and to a lesser degree also to changes in their 
amplitudes. Multivariate analyses show that these slow oscillatory phenomena covary with broad-band changes across multiple frequency bands 
(delta, theta, alpha, beta and gamma). Aside from such broad-band changes, there are also more specific interactions between slow and fast 
oscillations. Memory-related changes in gamma amplitude, for instance, covary with memory-related changes in theta amplitude and in theta phase 
alignment. Such theta-gamma covariance might be related to similar observations made in animal studies of information flow and coding in the 
medial temporal lobes. Finally, there are also prominent memory-related changes in theta amplitude that are topographically distinct from memory 
related changes in ERFs. Thus, aside from insights into the neural phenomena that underlie the memory-related changes in ERFs, the study of 
oscillatory rhythms in humans can reveal additional, more mechanistic insights into declarative memory and its links to animal studies. 
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EEG Oscillatory Phase, its Effects on Sensory Processing and its Modulation by Attention. 

Schroeder C.E., Lakatos P., Chen, C. & Karmos G. 

Cogn. Neurosci. Progr. Nathan Kline Institute, Orangeburg, USA; 

Dept. Neurosci., Albert Einstein Coll. Med. and Hungarian Acad, of Sci., Budapest 

Oscillatory activity ranging from very low frequency (delta) to very high-frequency (gamma) activity is ubiquitous in the brain, and its dynamic 
variations reflect changes in state that have important functional consequences. We examined the general question of an “oscillatory hierarchy” 
controlling cortical excitability, and its functional significance in terms of sensory processing and attentional modulation. Laminar current source 
density (CSD) and multiunit activity profiles were sampled from A1 with linear array multielectrodes, in awake macaque monkeys. Noteworthy 
qualities of spontaneous CSD include the periodic, low (theta) frequency CSD components in Layer 3, and higher (gamma) frequency CSD 
components in Layer 4. Instantaneous theta and gamma power and phase were extracted by applying the Morlet wavelet transformation to the CSD 
profile. Power was averaged for every 1° of phase to calculate the power-phase relation. We find that: 1) Delta, theta and gamma power are all largest 
in the supragranular layers. 2) Enhancement of gamma power and associated action potential bursts are “phase locked” to spontaneous, as well as 
stimulus-related theta oscillations. 3) Theta power is largest at a specific phase of delta oscillation. 4) Like gamma power, stimulus-evoked response 
amplitude depends on prestimulus theta phase. 5) Auditory stimuli cause partial phase re-setting of spontaneous delta (and thus theta and gamma) 
oscillations. 6) Delta phase-alignment facilitates interareal communication. 7) Though averaged auditory responses (ERPs) themselves reflect 
primarily stimulus-evoked neural activity (see, e.g., Shah et al., Cerebral Cortex, 14, 2004), attentional modulation of auditory ERPs appear to result 
mainly from partial phase resetting of ongoing neuronal oscillations. 8) Along with partial phase-resetting, auditory attention induces change in the 
frequency of theta band activity, so that it more closely approximates the average 1ST These effects have both parallels with and differences from, 
visual cortical effects within the same experimental paradigm. 
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Working memory: physiological characteristics of prefrontal information processing as revealed by 

BOLD fMRI 

J. H. Callicott 

Clinical Brain Disorders Branch, NIMH 

Augmenting basic imaging paradigms such as motor movement or visual processing, neuroimaging studies of short term or working 
memory (WM) have been a popular means of exploring individual and group differences in information processing strategies. While not the only 
area engaged by WM cognitive challenge, the dorsolateral prefrontal cortex (DLPFC) typically is engaged whenever such tasks demand some 
combination of maintenance, manipulation, updating, minimization of interference, or response selection based on the contents of WM. Even as 
fMRI methodology and cognitive challenges become more sophisticated, the central questions regarding the relationship between fMRI activation 
and the basic behavioral characteristics of working memory (capacity and efficiency) remain relevant to larger questions regarding the means by 
which humans utilize this key cognitive resource. These questions might be formulated as follows: 1. As WM performance falls in healthy subjects 
begin to make errors, what can one infer regarding capacity based on relative DLPFC activation? 2. When performance remains stable (i.e., within 
capacity), what can one infer regarding efficiency based on differences that arise in DLPFC activation? and 3. What inferences can one make 
regarding information processing strategies when comparing across groups that differ in either WM capacity or efficiency? While earlier BOLD 
fMRI work supplied rather straightforward answers to the first two questions (i.e., DLPFC activation decreases reflect reduced capacity while 
DLPFC activation increases reflect reduced efficiency), more recent work provides a more complex picture. Furthermore, quantitative labels such as 
"increased" or "decreased" activation are troublesome descriptors when WM paradigms based on experiments in healthy volunteers are used to 
generate hypotheses regarding clinical populations. Therefore, a review of the recent literature regarding prefrontal information processing as 
revealed by BOLD fMRI is a useful starting point as the field moves towards a greater appreciation of the temporal and spatial dynamics of WM 
offered by more sophisticated fMRI approaches combined with MEG. 


MEG activation comparison to fMRI BOLD for a working memory task 

R. Coppola 12 , J.H. Calicott 2 , T. Holroyd 1 , B.A. Verchinski 2 , S. Sust 2 , D.R. Weinberger 2 
'MEG Core Facility, NIMH, Bethesda, MD, USA. 2 Clinical Brain Disorders Branch, NIMH. 

Cognitive tasks have been explored with metabolic and electrophysiological activation measures. Here we examine the relation between BOLD fMRI 
and MEG to a working memory task in the same group of normal volunteers. Subjects performed a run of sixteen 22 second blocks alternating 
between a 2-back and a 0-back visually presented task. Stimuli were one of 4 numbers; response was a button press to match the current number (0- 
back) or the number from 2 trials previous. MEG was recorded on a CTF 275 channel system. Synthetic aperture magnetometry was used to analyze 
task related activation differences in theta, alpha, beta, and gamma frequency bands. SAM creates an optimum spatial filter from the covariance 
between the active state (2-back) and control state (0-back) to calculate a 3-d source image comparing the source strength for specified time windows 
for the two states. In this case we used time windows comprising either 500 msec after the stim, 500 msec centered on the response, or the middle 14 
sec of the block. Group t-maps were created in AFNI by aligning to Talairach the individual z-scored SAM volumes registered to their structural 
MRI. BOLD fMRI was collected at 3T (GE EPI-RT) to identical task design. Data were structurally aligned, smoothed, and normalized to MNI space 
using SPM99. Single suhj block design contrast maps for 2b>0b were followed by a group second level analysis and thresholded at approximately 
p<.001. This map showed activation of the working memory network as described in several studies. Of the 12 SAM maps (4 freq bands by 3 time 
windows) only the beta band for the response window showed very high agreement with the BOLD activation. At comparable p level threshold the 
MEG showed beta desynchronization with a very similar spatial extent to the BOLD activation. This relationship will be important for better 
understanding of the physiological basis of neuroimaging activation. 
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Working memory impairments after severe Traumatic Brain Injury: What does MEG reveal 

Asloun S (1-2)., Croize A.C.(2), Ducorps A.(3), Schwartz D.(3), Garnero L.(4), Azouvi P.(l) 

(1) : Neurological Rehabilitation unit, R. Poincare Hospital, Garches, France 

(2) : Inserm U.483, University of Paris VI, France 

(3) : MEG centre, la Salpetriere Hospital, Paris, France 

(4) : LENA-CNRS UPR 460, la Salpetriere Hospital, Paris; France 


ABSTRACT 

Many studies of working memory (WM), both in healthy subjects and severe Traumatic Brain Injury (sTBI) patients, with fMRI suggested the 
involvement of a prefrontal and parietal network activation. To our knowledge, no study of WM after sTBI has been reported with MEG imaging. 
This technique is however well suited to assess the slowing down of information processing observed in TBI patients. The aim is to assess the spatial 
and temporal profdes of brain activity during a WM task in sTBI compared to controls. Ten healthy subjects and five sTBI were submitted to a n- 
back task, with three conditions of increasing WM load. First results are in agreement with cortical networks described in previous fMRI studies. 
Both patients and controls showed a prefrontal and parietal network activation. The novelty of these results was to reveal that all stages of 
information processing, involving early visual sensory signals, were slowed down in the patients group. This slowing down increased with later 
cognitive stages of information processing. Furthermore, patients showed a higher signal intensity than controls. Finally, while controls showed an 
activation involving mainly the left hemisphere, activation was bilateral in sTBI. MEG appears thus to he an interesting technique for the study of 
impairments in WM after sTBI. It reveals that these patients showed both a global slowing down of information processing and a higher and more 
widespread level of activation, presumably as adaptative mechanisms. 

KEY WORDS 

Traumatic Brain Injury, working memory, n-back, MEG imaging. 

INTRODUCTION 

This study aimed to : assess the cerebral correlates of working memory impairments after severe TBI, explore the spatiotemporal brain dynamics 
using a working memory (WM) task and objectivize which stages of information processing are slowed. 

METHODS 

Five sTBI (mean age 25±3.69, educational level 15±3.20, score<8 on the Glasgow Coma Scale) and ten healthy subjects (mean age 24.5±3.92, 
educational level 15±3.07) were examined by MEG. 

A working memory (WM) task (n-back) was used in this study. A letters sequence were presented 0 Buck E l^tock 2 Duck 

on a screen during 500ms, with 3000ms interstimuli time. Subjects were trained before performing the ^ j q 

MEG imaging to make sure they understood well each condition of the task. There were three difficulty n R n X n K . 

levels of increasing WM load : 0, 1 and 2-Back (Figure I). In the 0-Back condition, subjects were ^ ■ 11 v X 

asked to press a “yes” key (keyboard) for a pre-detennined target letter (“S”). In the 1 -Back condition, n ' D C n ' 

they had to answer “yes” if the current letter was matched to the previous one. And for the 2-Back Figure 1 : Illustration of the three 

condition, they had to answer “yes” if the current letter was identical to that appeared 2 hack in the difficulty level of the n-back task 

sequence. 

For each participant, Event-Related magnetic Field (ERF) were recorded with a CTF (Vancouver, Canada) whole-head system (151 sensors). A 
3D T1 MRI acquisition (Philips, 1.5T) was also performed for superimposition for MEG acquisitions. For each subjects, 6 runs were recorded with 2 
conditions (0-1, 1-2, 0-2-Back) in the same run. Presentation order of runs was randomized over all subjects. Performances and Reaction Time (RT) 
of each subject were recorded. The signals of the brain were filtered with a 0.626-40 Hz bandpass. Eyes movements artefacts and cardiac artefacts 
were removed. Data were averaged according to each condition. ERF signals were analysed using a classical Minimum Norm Algorythms (MNE 
ERF) (Brainstrom 2001, S. Baillet, J.C. Mosher, R.M. Leahy). Z-scores calculations were performed over each MNE ERF with respect to a baseline- 
MNE signals (from -450 to -150ms). 

RESULTS 

The preliminary results were in agreement with cortical networks described in previous functional imaging studies. In the 2 groups, sTBI and 
controls, cortical networks found involved both frontal and parietal cortex. Both control subjects and sTBI showed (1) cortical activations in posterior 
parietal (BA39/40) and left prefrontal (BA 9/40) areas and (2) an increase of cortical activation intensity with the rise of WM load. However, for 
sTBI, signal intensity was higher than control subjects for each condition, especially in the prefrontal areas (bilateral). Furthermore, sTBI showed a 
disrupted activation pattern including supplementary areas (bilateral frontopolar areas, right temporal area). 

These results also showed a significant slowing of information processing in sTBI, compared to control subjects. This slowness involved the early 
visual sensory signal (115 ms for controls, 121.4 ms for sTBI) and increased with later cognitive processes. 

DISCUSSION 

Patients with severe TBI suffered from several cognitive deficits. Deficits in working memory and mental slowness are among the most frequent 
complaints. The novelty of this study was to assess this slowing down of information processing with MEG. Indeed, our results revealed that 
slowness involved every stages of information processing, including visual perception stages, but slowness seemed particularly marked for later 
stages involving the prefrontal cortex. 

Cerebral activation in both groups was found in similar regions of the frontal and parietal lobes. However, compared with controls, TBI patients 
showed a higher and more widespread signal intensity. This could be linked to differences in recruitment of resources or in solving strategy 
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differences between sTBI and controls subjects. It could also be related to a compensation of deficit activation in sTBI, that is, an increase of 
prefrontal activation could be necessary to maintain performance levels as high as control subjects. 

Previous studies in PET and fMRI showed a bilateral frontal and parietal activation. Severe TBI had supplementary areas activated, particularly 
the temporal areas. This could compensate a deficit in activation of cortical network needed to solve the n-back task by recruiting remote regions. 
The n-back task performances are known to depend on resources in working memory and temporal areas are known to be involved in long-term 
memory processes. The limitation of resources in working memory after sTBI could explain the involvement of long-term memory processes to 
perform the task. 

ACKNOWLEDGEMENTS 

This study is supported by a grant from Garches Institute, R. Poincare Hospital, Garches, France. 

REFERENCES 

Azouvi P, Jokic C, Van Der Linden M, Marlier N, Bussel B. Working memory and supervisory control after severe closed head injury. A study of dual 
task perfromance and random generation. Journal of Clinical and Experimental Neuropsychology, 1996, 18, 317-337. 

Fontaine A, Azouvi P, Remy P, Bussel B, Samson Y. Functional anatomy of neuropsychological deficits after severe traumatic brain injury. 
Neurology, 1999 Dec, 10, 53(9): 1963-8. 

Lewine JD, Davis JT, Sloan JH, Kodituwakku PW, Orrison WW Jr. Neuromagnetic assessment of pathophysiologic brain activity induced by minor 
head trauma. Am. J. Neuroradiol, 1999 May, 20(5):857-66. 

McAllister TW, Saykin AJ, Falshman LA, Sparling MB, Jonhson SC, Guerin SJ, Mamourian AC, Weaver JB, Yanofsky N. Brain activation during 
working memory 1 month after mild traumatic brain injury, Neurology, 1999, 53, 1300-1308. 

McDowell S, Whyte J, D'Esposito M. Working memory impairments after traumatic brain injury : evidence from dual task paradigm. 
Neuropsuchologia, vol 35, N° 10, pp 1341-1353, 1997. 

Christodoulou C, DcLuca J, Ricker JH, Madigan NK, Bly BM, Lange Q Kalnin AJ, Liu WC, Steffener J, Diamond BJ, Ni AC. Functional magnetic 
resonance imaging of working memory impairments after traumatic brain injury. J Neurol Neurosurg Psychiatry, 2001, 71:161-168. 


9 



S2 


Neuromagnetic Correlates of Memory Dysfunction: Lessons from the Clinic 

Jeffrey David Lewine, Ph.D.[l,3], Michael Funke, M.D., Ph.D. [1], Cathrine Chong, Ph.D. [1], Robert Thoma, Ph.D. [1,2], Erin 

Bigler, Ph.D., [1], Jose M. Canive, M.D.[2]. 

[1] Functional Brain Imaging Program, University of Utah School of Medicine; [2] Department of Psychiatry, Albuquerque Veterans 
Administration Medical Center; [3]Hoglund Brain Imaging Center, Departments of Neurology and Psychiatry & Behavioral Sciences, 

KUMC. 

Introduction: Memory dysfunction results in a significant compromise of quality-of-life for a wide range of patients, so there is growing interest in 
the development of better methods for early identification of mnemonic compromise and in the development of new therapeutic strategies 
(behavioral and phannacological). A major challenge in these endeavors is posed by the multifaceted nature of memory (event versus rule memory, 
implicit versus explicit, short versus long term) and the multistage nature of memory processes, different clinical populations showing different 
semiology or breakdown at different stages. Electrophysiological measures may be of particular utility because they may be more sensitive to the 
early breakdown of memory processes than behavioral measures. 

Materials and Methods: In a series of clinical studies, we have sought to identify both stimulus evoked and spontaneous activity correlates of 
mnemonic dysfunction in a range of patients including those with epilepsy, dementia, head trauma, and schizophrenia. 

Results'. In patients with memory problems subsequent to mild head trauma, abnormal low frequency activity is commonly seen over the anterior 
temporal lobes. Such slowing may be present, even when stimulus evoked signals appear to be intact. In patients with epilepsy and unilateral 
hippocampal sclerosis, abnormal spontaneous and abnormal evoked signals are seen, especially using tasks that tap into “event” memory, such as the 
running recognition task. Patients with schizophrenia, show evidence of impairment in both basic aspects of sensory evoked memory and working 
memory, often in the absence of low frequency spontaneous abnormalities. Finally, patients with dementia show abnormal spontaneous and evoked 
signals, especially at late stages of disease. 

Discussion and Conclusions: MEG is useful for assessing some of the neurobiological substrates of memory dysfunction in clinical populations. A 
multistage approach is useful for defining where breakdown occurs in a complicated information processing system. 
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From vibration in the ear to abstraction in the head: early auditory analysis 

D.Poeppel 

University of Maryland College Park, USA 

A central question for research on speech perception concerns how to transfonn input signals that are continuous (wavefonns at the auditory 
periphery) into the discrete elements that underlie lexical representation (distinctive features, segments, syllables). The information arriving at the 
periphery must be fractionated in to chunks of the appropriate granularity, and the necessary features must be extracted to permit word recognition 
and the recognition of prosody. Acoustic features that play a central role for the construction of speech sounds are frequency, bandwidth, and 
temporal stimulus structure. Recent MEG studies are discussed that exemplify how these features are being studied. In particular, we consider results 
(i) on frequency representation in space (tonotopy) and time (tonochrony), (ii) on the sensitivity to noise and bandwidth, and (iii) on temporal 
integration of auditory information. The recognition of complex sounds and the distinction between speech and non-speech must ultimately build on 
these mechanisms. 

To illustrate one domain in more detail, results on the sensitivity to noise and bandwidth modulation are described. Complex sounds, specifically 
natural speech sounds, are broadband signals that have internal structure (e.g. formants) with complex spectro-temporal dynamics; therefore a better 
understanding of the processing of simple signals with bandwidth is necessary. MEG results demonstrate a high degree of sensitivity to bandwidth 
modulation in early cortical responses. For example, increasing bandwidth systematically alters the classical M50 and M100 evoked responses and 
also systematically decreases the steady-state responses. Such results identify physiological properties that form the basis for constructing auditory 
objects that become speech sounds and ultimately interface with stored lexical representations. 

DP is supported by NIH R01 DC05660. 


Perception and representation of speech-specific information 

C. Eulitz 

RWTH Aachen and University of Konstanz, Germany 

No word is ever pronounced alike twice. The resulting amount of acoustic/phonetic details in the sensory representations of speech has to be reduced 
to an optimal set of manageable information by transforming them into phonological categories. These categories are abstract and acoustic 
differences among category members are irrelevant. This makes the higher level representations of speech sounds sparse and good for computation. 
Recent MEG studies are discussed that exemplify how phonological features (the smallest building blocks of language) are extracted from speech. It 
will be demonstrated that the NIOOm indexes the activation of cortical structures tuned for the extraction of phonological features. Spatially distinct 
cortical structures are involved in the extraction of the mutually exclusive phonological features and show moreover a differential timing of 
activation. 

The fine structure of the higher level representations of speech sounds can be studied by means of the mismatch negativity/field. Recent studies are 
discussed that demonstrated how the mismatch response can be used to learn more about the functional organization of the mental lexicon. 

Research was supported by the Deutsche Forschungsgemeinschaft. 
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Sensory processing in dyslexia 

H. Renvall 

Brain Research Unit, Low Temperature Laboratory, Helsinki University of Technology, Finland 

Dyslexic individuals encounter difficulties that often extend beyond the skills directly needed for reading. For example, language-leaming-impaired 
children are slow in processing sounds presented in rapid succession or containing fast frequency transitions. These auditory deficits, at a time scale 
of up to a few hundreds of milliseconds, persist to adult age. Recent behavioral studies have implied that dyslexic subjects have temporal processing 
deficits in other sensory modalities as well. 

We have applied magnetoencephalography (MEG) to characterize the temporal impairment in auditory and tactile domains in dyslexic adults. Our 
results demonstrate that dyslexic adults are deficient in processing sounds and acoustic changes presented in rapid succession within tens to hundreds 
of milliseconds. In line with the proposed pansensory processing deficit, responses to rapidly presented tactile stimuli were diminished in the right 
somatosensory cortex of dyslexic subjects. On the basis of these and earlier findings we have proposed that limitations in both modality-specific and 
more global attentional capacities could prolong input chunks for all senses, thereby leading to anomalous cortical representations in dyslexic 
subjects [1], MEG provides a powerful tool for characterizing sensory processing deficits in dyslexia, but further studies are needed to clarify the 
relationship between sensory-specific cortical processing and phonological processing deficits in dyslexic subjects. 

[1] Hari, R„ Renvall, H. 2001. Impaired processing of rapid stimulus sequences in dyslexia. TICS 5, 525-532. 


Early cortical processes in fluent and impaired reading 

R. Salmelin 

Brain Research Unit, Low Temperature Laboratory, Helsinki University of Technology, Finland 

MEG studies have demonstrated a sequence of two distinct processes in word perception within the first 200 ms. Low-level visual feature analysis 
occurs at 100 ms, around the occipital midline, and it is followed by letter-string-specific (LS) activation at 150 ms, in the left inferior 
occipitotemporal cortex [1], The source area of the LS response is located slightly posterior (BA19/37) to the Visual Word Form Area (VWFA) 
identified in hemodynamic studies (BA 37) [2,3], The LS response does not differentiate between words and nonwords but VWFA apparently does. It 
is plausible that the early transient MEG signal reflects the onset of letter-string-specific analysis which is not detected in, or does not dominate, the 
fMRI/PET signal. More anterior areas along the ventral stream may then be recruited where neurons are increasingly sensitive to the word-likeness of 
the letter-strings but show weaker synchronization or are less rigorously time-locked to stimulus presentation and would thus not be detected with 
usual MEG analysis techniques [4,5], In dyslexic individuals, the non-stimulus-specific visual feature analysis is intact but the reading process is 
disrupted at the stage of letter-string-specific analysis, which is non-detectable or abnormally weak [6]. Face-specific analysis in the same area and 
time window, however, appears to be normal [7]. Underactivation of the left inferior occipitotemporal cortex in dyslexia is the most consistent 
finding in hemodynamic recordings, as well [8], 

[1] Tarkiainen et al. 1999. Brain 122, 2119-31. [2] Tarkiainen et al. 2002. Brain 125, 1125-36. [3] McCandliss et al. 2003. Trends Cogn Sci 7, 293-9. 
[4] Cornelissen et al. 2003. J Cogn Neurosci 15, 731-46. [5] Salmelin et al. 2004. Sci Studies Reading 8, 257-72. [6] Helenius et al. 1999. Cereb 
Cortex 9, 476-83. [7] Tarkiainen et al. 2003. Neuroimage 19, 1194-204. [8] Paulesuet al. 2001. Science 291, 2165-7. 
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Early lexical and syntactic processing as revealed by the Mismatch Negativity: IMMNm and sMMNm 

Friedemann Pulvermuller & Yury Shtyrov 
MRC Cognition and Brain Sciences Unit, Cambridge 

If a syllable completes a meaningful word, it elicits a larger Mismatch Negativity (MMN) than in a context where it is part of meaningless 
pseudowords. The main source of the magnetic MMN to lexical items localized in left superior temporal cortex became active at —150 ms after 
stimulus information allowed for unique word recognition. The superior temporal activation was followed, with a delay of ~20 ms, by left inferior 
frontal activation, indicating that lexical processing is reflected by an ordered spatio-temporal pattern of activity, the near simultaneous ignition of at 
least two distinct left perisylvian sources. Grammaticality of a word sequence was also reflected by dynamics in the MMN. Words presented in a 
context where they violate agreement rules led to larger magnetic MMNs than the same words in grammatical context. The time ranges of the early 
lexical and syntactic dynamics of the MMN were the same, supporting early multilevel access to linguistic information. Furthermore, our results 
show that early lexical and syntactic processes do not require focussed attention. In all experiments, acoustic properties of the stimuli were exactly 
controlled and we hypothesize that this is an important feature necessary for revealing the early (200 ms and below) neurophysiological 
manifestations of lexical access and syntactic processing. 

[1] F. Pulvermuller el al., Neuroimage 14, 607-616 (2001); Y. Shtyrov, F. Pulvermuller, NeuroReport 13, 521-525 (2002); F. Pulvermuller, Y. Shtyrov, 
Neuroimage 20, 159-172 (2003). 

[2] F. Pulvermuller, Y. Shtyrov, R. J. Ilmoniemi, Neuroimage 20, 1020-1025 (2003); Y. Shtyrov, F. Pulvermuller, R. Naatanen, R. J. Ilmoniemi, 
Journal of Cognitive Neuroscience 15, 1195-1206 (2003). 

[3] R. Assadollahi, F. Pulvermuller, NeuroReport 14, 1183-7 (2003); O. Hauk, F. Pulvermuller, Clinical Neurophysiology 115, 1090-103 (2004). 


Lexical Access and Beyond: The M350 

A. Marantz 

Department of Linguistics & Philosophy, MIT, USA 

Recent MEG work reveals that visual word presentation evokes a response from a temporal source near auditory cortex peaking around 350ms that 
indexes “lexical access” as understood in the cognitive literature [1], I review evidence from the KIT/MIT MEG lab and others supporting the M350 
as the first response component reflecting access to the representations of roots of lexical items connecting sound and meaning (e.g. also [2]). 
However, apparent word frequency effects at around ~150ms from occipito-temporal “visual letter-string” areas point to the presence of lexical 
representations earlier in word recognition and in specialized visual processing areas (e.g., [3]). Evidence will be presented from masked priming 
experiments that the MEG Ml70 letter-string response is sensitive to the morphological structure of a word (performs affix-stripping) but not 
necessarily to the lexical status of a letter string (pseudo-affixes are stripped, e.g., the -er of brother). The evidence supports the hypothesis that the 
specialized letter-string brain region is sensitive both to orthographic frequencies (and letter-sequence frequencies) and to morphological structure, 
but not directly to lexical properties of word root representations. I will emphasize the convergence of MEG, ERP, fMRI and behavioral evidence on 
the time-course of visual word processing and on the representations accessed by different brain regions. 

[1] Pylkkanen, L., Marantz, A. 2004. Tracking the time course of word recognition with MEG. TiCS 7.5, 187-189. 

[2] Friedrich, C.K. et al. 2004. ERPs Reflect Lexical Identification in Word Fragment Priming. JCNC 16.4, 541-552. 

[3[ Sereno, S.C., Brewer, C.C., O’Donnell, P.J. 2003. Context Effects in Word Recognition. Psychological Science 14.4, 328-333. 
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Semantic Processing of Words and Sentences 

Ksenija Marinkovic 

Martinos Center for Biomedical Imaging, Mass. General Hospital, Harvard Med. School, Boston, USA 

Understanding language relies on concurrent activation of multiple areas within a distributed neural network. In agreement with other neuroimaging 
evidence, studies using the anatomically-constrained MEG approach suggest that both written and spoken words are initially processed in their 
respective ventral modality-specific streams. The subsequent stage reflecting semantic access and integration (“N400”) is subserved by distributed 
brain regions that are simultaneously active for a protracted period of time, culminating in a generic process of word comprehension. Predominantly 
left temporal, prefrontal and medial prefrontal areas may provide specialized contributions in an interactive and mutually dependent manner during 
the process of constructing the meaning that fits best in the given context. Different lines of evidence suggest that lexical access, semantic 
associations, and contextual integration may be simultaneous and interactive as the brain uses available information in a concurrent manner, with the 
goal of rapidly comprehending verbal input. The right prefrontal cortex additionally participates in certain aspects of contextual integration such as 
comprehending words with weak semantic associations, alternative meanings, or jokes. Indeed, jokes seem to engage the right prefrontal cortex 
following the N400, during the phase of retrieving the alternate meaning so that the “twist” can be incorporated into the joke context. Because the 
same areas may participate in multiple stages of verbal processing, it is crucial to consider both spatial and temporal aspects of their interactions to 
appreciate how the brain derives meaning from words. 

Supported by NIH (AA13402, P41RR14075), and the MIND Institute. 


Representation of meaning in the brain 

B. Rockstroh, R. Assadollahi, A.Low, T. Elbert 
Department of Psychology, University of Konstanz, Germany 

We used Magnetic Source Imaging to elucidate the representation of different semantic and syntactic elements in a series of studies - all exploiting 
event-related responses to visual stimuli that induced or were associated with syntactic or semantic processes. Representation of meaning was 
examined through pictures (Exp. 1) or words (Exp.2) of objects from superordinate semantic categories (e.g. flowers, animals). In response to 
pictures, the spatio-temporal correlation of brain activity was larger for stimuli within than across superordinate categories around 200 ms in the 
right-parietal-temporal lobe, and around 400 ms in the left temporal lobe. In response to words proximity of representations within a category was 
indicated by a smaller Euclidean distance between dipoles, relative to dipole distances across categories; this focal topography distinguishing 
categories was prominent around 200 ms post-stimulus and only within the left temporal lobe. Results indicate that stimuli are mapped onto semantic 
categories along a hierarchical stream with its endpoint in the left temporal region. The interface between syntax and meaning was examined by 
varying the verb’s argument structure in a series of 390 verbs requiring one (e.g. Peter , snored’.), two (e.g. Peter saw John.) or three (Peter gave John 
s.th.) obligatory arguments. Left-temporal activity at 250 ms distinguished the number of arguments, whereas processing the same verbs within a 
context (e.g. Peter snored.) activated areas in the left inferior frontal gyrus, indicating the distinction of the verb’s representation from the processing 
of its argument structure. 

In conjunction with the literature, we conclude, that syntax and meaning are represented in macroscopically different areas (networks) within the left 
hemisphere. Feature selection along the level of abstractness may occur through top-down influences on hierarchical streams of information 
processing. While semantics reaches its categorical stage within the temporal lobes, the highest stage of syntax analyses may be reached in the left 
inferior frontal region (possibly Broca’s area). 

Research was supported by the DFG and the Volkswagen Stiftung. 
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Clinical Multichannel MCG in Unshielded Hospital Environment 

Brisinda D., Meloni A.M., Fenici R.. 

Clinical Physiology -Biomagnetism Center, Catholic University, Rome - Italy 


ABSTRACT 

From November 5 th , 2001 to May 19 th , 2004, 545 patients (177 with arrhythmias, 67 with WPW syndrome, 60 with Ischemic Heart Disease 
(IHD), 129 with different kinds of cardiomyopathy, 106 normals, 6 FMCG) have been consecutively investigated at the Catholic University of Rome, 
with unshielded Multichannel Magnetocardiographic Mapping (MMCG): 20 with the 9-channel system only and 525 with the 36-channel system 
(207 of them with both systems). 107 patients were investigated also after physical stress, carried out with a standard bicycle ergometer. In all 
patients MMCG was recorded at least three times, to check for reproducibility and/or for clinical follow-up, for a total of more than 1600 recordings. 
METHOD: MMCG was performed, with both the 9-channel and the 36-channel systems, at 1 kHz in the bandwidth DC-100 Hz. In the last 200 pts, 
12-lead ECG was simultaneously recorded with amagnetic electrodes. On each patient file, post-processing and signal analysis for the quantitative 
assessment of ventricular repolarization and for 3D localization and electroanatomical imaging of cardiac arrhythmias, were carried out 
independently with two different approaches and software programs developed by CMI and by Neuromag (Finland). RESULTS: The results with the 
two methods have been compared. For 3D electroanatomical integration of MMCG localization results, 3D cardiac models have been used, 
constructed from patient MRI and/or from orthogonal fluoroscopic images taken at the moment of MCG recording. CONCLUSIONS: Qualitative 
reproducibility of MMCG was satisfactory. However the estimate of quantitative parameters has shown a certain degree of variability, which deserves 
further evaluation. 

KEY WORDS 

Multichannel Magnetocardiography, Body Surface Mapping, Ischemic Heart Disease, Arrhythmias, WPW Syndrome, Cardiomyopathy, Fetal 
Magnetocardiography, Stress Test. 

INTRODUCTION 

Previous research has shown that, compared to standard ECQ multichannel magnetocardiographic mapping (MMCG) provides non-invasive 
evaluation of cardiac electrogenesis, with similar investigation time, but higher spatial and temporal resolution [Tavarozzi, 2002], The diagnostic 
potentiality of MMCG spreads from three-dimensional (3D) electroanatomical functional imaging and localization of focal or re-entry 
arrhythmogenic mechanisms [Fenici, 1998] [Moshage, 1996] [Giorgi, 2004], with possible correlation of the site of origin of the arrhythmia with the 
abnormality of the anatomic substrate [Leder U, 1998], to quantitative estimation of ventricular repolarization (VR) abnormalities in patients with 
ischemic heart disease (IHD) [Hanninen, 2000] [Tsukada, 2000] or with other cardiomyopathy (CMP) [Shiono, 2003] [Agrawal, 2001], 

However, in spite of such interesting features, the high costs of the instrumentation and the prevailing opinion that MMCG could be performed 
only in magnetically shielded rooms have prevented the acceptance of this method as a clinical tool. More recently, after the development of 
multichannel systems, easy to install and reliable in any unshielded hospital (even emergency) environment at relatively low cost, the number of 
clinicians interested in MMCG is rapidly increasing. MMCG is now under investigation as useful alternative for rapid screening, at rest, of patients 
with suspected IHD and for the “triage” of emergency patients with acute chest pain and still negative or dubious ECG and cardiac enzyme patterns 
[Steinberg, 2004] [Tolstrup, 2004], 

All studies published so far have been carried out on relatively small groups of patients, and with different methods for data reduction and 
analysis. Therefore the possibility to compare and/or pool data coming from different centers is limited. This might be particularly evident and 
important for large-scale applications such as for the detection of cardiac ischemia. Indeed the analytic methods and parameters used for the 
diagnosis of ischemia with MCG are not univocal and do surely need a better definition and validation. Unfortunately the mechanisms underlying VR 
abnormalities detected with MMCG in IHD patients are not fully understood yet [Tsukada, 2000], and it is well known that similar VR abnormalities 
have been observed in patients affected by different kind of CMP (hypertrophic, dilative or inflammatory). 

At the Catholic University of Rome, after preliminary testing carried out with a 9-channel system [Fenici, 2001], we succeeded in installing a 
multichannel system for MMCG which provides good quality 36-point recording in only 90 seconds, even working in the unshielded catheterization 
laboratory [Fenici, 2003], 

In this paper we have reported our experience, working with this unique installation from November 5 th , 2001 to May 19 th , 2004. 

METHODS 

545 patients (pts), ranging in age from 38 months to 92 years (177 with arrhythmias, 67 with WPW syndrome, 60 with IHD, 129 with different 
kinds of cardiomyopathy, 106 normals), were consecutively investigated from November 5 th , 2001 to May 19 th , 2004, in the unshielded room 
equipped for invasive cardiac diagnostic and interventional procedures. To test the feasibility of unshielded fetal magnetocardiographic mapping 
(FMCG), 6 normal pregnant women were investigated, at different gestational epochs, from the 32 nd gestational week. In all patients MCG was 
recorded at least three times, to check for reproducibility and/or for clinical follow-up, for a total of more than 1600 recordings. 

MMCG was initially recorded with the “scout” 9-channel system (227 pts) and then with the 36-channel system (525 pts; 207 of them previously 
studied with the 9-channel) (CardioMag Imaging Inc., Schenectady, NY) (CMI). The latter system features 36 LT-DC-SQUID (Superconducting 
Quantum Interference Device) sensors coupled to second order gradiometers (baseline: 50-70 mm) with pick-up coils diameter of 19 mm and sensor- 
to-sensor spacing of 40 mm. The distance between the measuring sensors, kept at liquid helium temperature and arranged on the same plane, and the 
flat bottom surface of the cryostat is 19 mm. With the automatic electronic noise suppression system (ENSS), the instrumentation features a 
sensitivity of about 20 fT/Hz!4 at 1 Hz, with balance stability of gradiometers better than 0.01%. MCG signals and one reference ECG lead are 
simultaneously recorded at the sampling rate of 1 kHz, in the bandwidth from DC to 100 Hz. In the last 200 pts, 12-lead ECG was recorded 
simultaneously using amagnetic electrodes and cables ( GE-Prucka CARDIOLAB System), 107 of them (with IHD, essential hypertension, or normal 
controls) were investigated also before and 1,3,5 and 7 minutes after physical stress with a standard bicycle ergometer. 


15 



S 5 


For the quantitative assessment of MCG parameters and for 3D localization and electroanatomical imaging of cardiac arrhythmias, each patient 
file, is analyzed independently with two different software packages, developed by CMI and by Neuromag (Finland), respectively. For 3D source 
localization the Effective Magnetic Dipole (EMD), the Equivalent Current Dipole (ECD), Current reconstruction (CR) and current density imaging 
(CDI) models, are available. 

For 3D electroanatomical integration, MCG localization results are transferred into the patient’s MRI, and/or 3D elastic models of the heart 
adapted to patient’s geometry inferred from orthogonal fluoroscopic images taken at the moment of MCG recording. 

To assess VR, MCG parameters proposed by different authors [Hanninen, 2000] [Steinberg, 2004] [Tolstrup, 2004] [Brisinda, 2003] are measured 
and compared: 1) magnetic field (MF) orientation (STa angle and Ta angle) [Flanninen, 2000]; 2) three MF dynamics parameters [i.e. +/- poles: 
angle (A), distance (D) and ratio (R)] during the during the T-wave interval; 3) quantitative scores of the T-wave (EXT, Machine learning (ML) and Q 
scores) were also automatically calculated with the "Machine Learning” approach [Embrechts, 2003], 

The same patient file were analyzed after preprocessing with different filtering modalities (such as digital low-pass filtering at 20 Flz, or digital 
adaptive filtering of the 50 Hz power line noise only) [Fenici, 2004], 

To improve the S/N ratio of FMCG, Fixed Point Independent Component Analysis algorithm (Fast-ICA) was to separate the maternal and fetal 
signals in order to reconstruct the time course of the electrical activity of the fetal heart alone from raw FMCG data. ICA permits a completely blind 
separation of signal sources starting from given mixed input signals under the condition that the mixed traces have components that are independent 
from one another [Hyvarinen, 1999], In two fetuses full mapping of cardiac magnetic field was possible. 
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RESULTS 

As for the 9-channel system [Fenici, 2001], the “radiofrequency” 
interferences have represented, the major problem during the 36- 
channel installation. This problem is nowadays fully overcome and 
the system works reliably 24 hours a day, with identical sensitivity 
even during rush hours. 

The preliminary results of clinical applications of the 36-channel 
system have been described in different papers. Here we will only 
briefly summarize few more relevant results. 

For instance MCG imaging and localization is particularly useful 
to properly classify accessory pathways, especially those localized in 
septal area, which determine complex activation patterns during the 
delta wave, of difficult interpretation with ECG only. Our results, 
obtained comparing MCG and basal ECG, have showed that MCG 
localization of ventricular preexcitation was more precise (92.8% vs 
52,3%) and allowed to identify dual accessory pathway not visible at 
ECG [Brisinda, 2004a], 

MMCG sensitivity was good enough to detect VR abnormalities 
in patients with chronic IHD, using the same analytic method 
previously reported by the Helsinki authors [Hanninen, 2000], both a 
rest, in which it evidenced morphological abnormalities of the MFD 
during the ST-segment in 88% of pts, and also during effort test, with 
a rotation of the ST a angle toward abnormal values in all IHD 
patients, in spite of normal ECG (Figure 1). 

The sensitivity of the MF dynamics parameters A, D, and R in 
detecting IHD was 32%, 42% and 42% with 20 Hz LPF, and 47%, 
74% and 63% with 50 Hz AF respectively. The difference was 
statistically significant for parameter D only (p<0.01). The sensitivity of the three automatic scores calculated with the CMI software (EXT, ML and 
Q, with 20 Hz LP) was 63%, 74% and 79% respectively [Fenici, 2004], 

MCG has good enough sensitivity to detect electrophysiological abnormalities induced by essential hypertension. As compared to controls three 
major differences were found: 1) prolongation of repolarization intervals, 2) increment of QT dispersion 3) abnormal orientation (a angle) of the ST 
MF. 
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Figure 1. Example of abnormal MF distribution and dynamics during 
the T-wave, in spite of nonnal ECG at rest (a). At the first minute 
after the termination of exercise (Rl’) the rotation of T-wave MF is 
evident (b), although VR at the ECG is still within nonnal limits. 


High-resolution analysis of atrial MF is feasible with unshielded MMCQ providing clear-cut differentiation between atrial depolarization and 
repolarization phenomena occurring during the P wave, which is traditionally considered due to atrial depolarization only. 

In agreement with previous human ECG studies [Bidoggia, 2000], gender-related differences of VR parameters were found also with MMCG. In 
particular, measurements of duration of total repolarization and of transmural dispersion of repolarization (through the measurement of Q-T peak , Q- 
Tend, J-Tpeak, J-T en d> T peak -T end intervals) were significantly longer in females. 

FMCG was possible only after the 30 lh gestational week. However the quality of fetal traces, recorded after the 32 nd gestational week, 
reconstructed with FastICA algorithms [Hyvarinen, 1999], was close to that of fetal signals obtained in the shielded room [Van Leeuwen, 2000], and 
good enough to identify the P, QRS, and T waves, to measure cardiac intervals and evaluate heart rate variability parameters both in the time and 
frequency domains [Brisinda, 2004b] (Figure 2). 

Finally, in spite of the unfavorable environmental condition contactless MMCG was feasible and reproducible with both the 9- and 36-channel 
systems, also for study of small animals, such as rats and hamsters, with minimum animal weight about 250 grams, for detection of both atrial and 
ventricular MF [Brisinda, 2004c], 
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DISCUSSION 

Although MMCG in shielded room remains a goldstandard for 
quality and experimental study of MCG recordings, it is nowadays 
evident that the trend is toward the use of unshielded instrumentation, 
easy to install and reliable even in emergency room. This new possibility 
is capturing the interest of clinicians worldwide with new installation in 
Unites States, Europe and China. In this study we have summarized our 
present experience with MMCQ in our unshielded hospital setting, where 
since the late seventies we tried to bring magnetocardiography at the 
patient’s bedside [Fenici, 1980], At that time one MCG mapping required 
more than two hours and pre-processing and data analysis took at least 
three additional hours. Our 36-channel system provides real-time 
monitoring of the cardiac magnetic field distribution during mapping and 
automatic data analysis in less than five minutes including source 
localization and VR assessment. 

Signal averaged MCG provided a clear-cut differentiation between 
atrial depolarization and repolarization, during the P wave. Interestingly 
such MCG parameters might be independent predictors of recurrences of 
paroxysmal atrial arrhythmias [Giorgi, 2004], The accuracy of MCG 
localization of paraseptal accessory pathways was validated with 
transcatheter ablation [Brisinda 2004a] and, in some patients, post¬ 
ablation mapping was also useful to understand the reason for 
interventional failure [Brisinda 2004d]. 

The quality of FMCG recorded in our unshielded hospital laboratory 
allowed reliable data analysis according to the recommended standards 
[Grimm, 2003] only after the 32 nd gestational week. This might obviously 
limit the clinical application of unshielded FMCG to the last two months 
of gestation, however provides the possibility to monitor fetal heart rate and well being on a simple outpatient basis, without any discomfort for the 
mothers. 

In patients with chronic IHD, the predictive value of rest MMCG was in the order of 70 to 80% as a function of inclusion criteria. This result is in 
good agreement with recent findings of other authors working with unshielded MCG [Steinberg, 2004] [Tolstrup, 2004] [Embrechts, 2003]. 
Furthermore we demonstrated that exercise-MCG is feasible with a conventional ergometer, as a complement to standard 12-lead ECG stress test. 
Stress-MCG improves the sensitivity of method to detect MF alterations due to exercise-induced myocardial ischemia. We could reproduce the 
typical MF rotation observed in previous studies carried out in magnetic shielding rooms [Hanninen, 2000] (figure 1). We did not test so far other 
parameters, such as the QRS integral, which have been recently reported to be sensitive to exercise-induced ischemia [Kanzaki, 2003], 

To conclude, as concerns the use of MMCG for the early (or emergency) detection of cardiac ischemia, our opinion is that, although the 
sensitivity of MMCG in detection VR abnormalities might better than that of rest ECG, further work is needed to define new analytic criteria to 
differentiate VR abnonnalities due to myocardial ischemia from those due to other cardiomyopathy, such as myocarditis, left ventricular hypertrophy, 
or dilatation [Shiono, 2003], In order to reach this goal experimental studies with animal models will be necessary. In this view we have successfully 
investigated the feasibility of MCG mapping of small animals with the same MMCG system used for clinical recording [Brisinda, 2004c], 

Finally, a further significant step forward for clinical application of magnetocardiography would be the construction of high-sensitive 
instrumentation based on second-order gradiometers coupled to low-cost non-cryogenic optical magnetometers [Bison, 2003]. 
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Figure 2. Example of unshielded fetal MCG mapping, (a) Real 
time recordings from top to bottom: maternal ECG FMCG row 
data, unfiltered FMCG after Fast ICA, smoothed FMCG, (b) 
Averaged FMCG, (c) Averaged FMCG waveforms recorded from 
36 points, (d) Fetal MF reconstruction of the peak of the QRS and 
source localization with the EMD model. 
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ABSTRACT 

The noninvasive diagnosis of cardiac ischemia related to coronary artery disease (CAD) remains a clinical challenge. Despite the large number of 
methods available in clinical routine, the predictive value of each of these techniques is still limited. Magnetocardiography (MCG) as a completely 
noninvasive approach offers new insights in the electrogenesis of the disease. In the course of the last decade a number of studies using 
biomagnetometers in shielded or unshielded environments have dealt with the identification of CAD patients diagnosed as such using coronary 
angiography as gold standard. As the availability of these systems in clinical use is still limited at present, studies have focused either on exercise- 
induced ischemia or on chronic ischemia and the infarct scar at rest. 

Different parameters have been developed based primarily on signal morphology, time intervals, source parameters or magnetic field map analysis. 
With respect to signal morphology, main work concentrates on ST-depression, ST-T signal amplitude as well as QRS and ST-T integrals. Dealing 
with time intervals most studies focus on QT interval. The evaluation of spatial distribution of QT dispersion in MCG, reflecting regional 
heterogeneity of repolarization, improved the identification of CAD patients at rest and during stress. Beside the calculation of the equivalent current 
dipole during de- and repolarization, parameters of the magnetic field orientation were used to identify CAD patients and localize exercise-induced 
ischemic regions. Heart rate adjusted alteration in the magnetic field orientation allowed the quantification of ischemia-induced changes in MCG. 
The estimation of current density (CDV) further enabled to separate healthy subjects from CAD patients at rest. In the course of interventional 
therapy CDV map distribution returned toward that of healthy subjects. In conclusion, several promising approaches exist to justify the application of 
MCG in the diagnosis of CAD for routine clinical use. 

KEY WORDS 

coronary artery disease, magnetocardiography, signal morphology, time intervals, source parameters, magnetic field maps, shielding 

INTRODUCTION 

Coronary artery disease (CAD) is still the leading cause of death in adults in industrial countries with a high annual incidence of the disease. For 
this reason scientific work emphasizes the early identification of these patients. Although coronary angiography is the gold standard at present, it 
should be restricted to those patients in whom interventional therapy appears to be indicated in order to avoid x-ray exposure and possible serious 
side effects. Despite the development of novel noninvasive diagnostic tests including new imaging techniques on the basis of magnetic resonance 
imaging or computed tomography, 12 lead ECG echocardiography and radionuclide imaging at rest and under stress are still the most common 
methods in routine clinical use. Despite technical advances, the predictive value of each of these approaches is still limited. Magnetocardiography 
(MCG), a completely non-invasive method without any side effects, allows the recording at the body surface of the magnetic fields generated by the 
electrical activity of the heart. MCG recording systems permit spatially and temporally accurate measurements of the very weak magnetic fields 
produced by currents flowing within myocardial fibers during cardiac activity. Beside the higher sensitivity of MCG to tangential currents in the 
heart compared to ECG it is also sensitive to vortex currents, which cannot be detected by the ECG. For these reasons, MCG can be regarded as an 
alternative non-invasive approach to cardiac electrophysiological phenomena, complementary to the electrocardiographic approach with a possibly 
improved diagnostic potential. 

The clinical value of MCG in the context of CAD covers several aspects: screening for subjects with high risk for the disease, the early diagnosis 
of patients with hemodynamically relevant stenoses of the coronary arteries with and without myocardial infarction (MI), the follow up of patients 
after interventional therapy and, lastly, the identification of CAD patients at risk for malignant arrhythmias. The following survey will focus primarily 
on the value of MCG in the diagnosis of the disease with respect to critical stenoses of the coronary vessel as documented by angiography. Various 
parameters derived from the magnetic signal and reconstructed magnetic fields under stress or resting conditions will be considered (Table 1). 
Furthermore the effect of shielding on the MCG procedure will be shortly discussed. 

ANALYSIS OF THE MCG SIGNAL 

Myocardial ischemia induces changes in the electrophysiological properties of the myocardium, resulting in a decrease in resting membrane 
potential and conduction velocity, with a dispersion of the activation wavefront. Transient acute ischemia causes well-recognizable changes in the 
ECG on the ST segment and the T wave. Performing the first single-channel MCG in a CAD patient after exercise testing, Saarinen et al. could show 
a depression of the ST-segment in the MCG signal. The ratio of the ischemic ST segment depression to R-wave amplitude after stress was even 
greater in the MCG than in the ECG [Saarinen, 1974]. Cohen et al. found TQ baseline elevation and ST depression after a two-step exercise test of a 
patient with CAD [Cohen, 1983]. Hanninen et al. could show a decrease of ST amplitude, ST slope and T-wave amplitude under exercise-testing in 
44 patients with ischemia documented by coronary angiography and exercise thallium scintigraphy. The optimal sites for the measurement of these 
parameters were over the abdomen. The reciprocal increase was found over the left parasternal area [Hanninen, 2002], 

DETERMINATION OF TIME INTERVALS 

QT interval prolongation has been associated with increased risk of malignant arrhythmias in patients after MI [Schwartz, 1978], The predictive 
value could be improved by the calculation of QT dispersion (QTd) measured as the difference between the maximum and minimum QT duration in 
the 12 lead ECG. QTd is assumed to reflect inhomogeneity of ventricular repolarization in different regions of the myocardium. These malignant 
arrhythmias are often ischemic in origin indicating that the increase in QTd could be due to myocardial ischemia [Higham, 1995]. Ideally, the 
measurement of QT intervals at numerous locations over the thorax - feasible with multichannel MCG - should permit a more detailed examination 
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of the spatial aspects of QT duration disparity. Indeed a better separation between CAD patients and healthy subjects was given on the basis of 36 
channel magnetocardiography compared to 12 lead ECG when the temporal approach - as expressed by QTd - was applied. Consideration of the 
spatial distribution of QTd improved its sensitivity in the detection of CAD. Increased QT variability between neighbouring sites and an alteration in 
the global repolarization pattern compared to healthy subjects could be displayed in CAD patients not only under stress [Hailer, 1999] but also in a 
resting condition [Hailer, 1999b, van Leeuwen, 2003], For quantification of these spatial differences, two smoothness indexes were calculated: SI, 
which quantifies the mean dispersion at each registration site and Sin, which normalizes SI by taking the overall spatial dispersion pattern into 
account (figure 1) [van Leeuwen, 1996]. Expanding the coverage area by using a 61 channel magnetometer that covers most of the thorax, the 
identification of CAD patients confirmed previous results [Klein, 2002], 


SOURCE PARAMETERS 

The ability to study the three 
dimensional localization of cardiac 
electrical activity makes MCG a potentially 
useful tool in the evaluation of the presence 
and localization of CAD. In the field of 
chronic ischemia and myocardial viability, 
current densities for the endocardial surface 
of the left ventricle were calculated and 
applied to predefined myocardial geometry 
acquired from MR images. The areas 
corresponding to the infarcted segments of 
the heart defined by the clinical reference 
methods revealed markedly decreased 
regional current densities [Leder, 1998], 
Current density estimation (CDE) was also 
tested in the localization of ischemia in 
patients with single and three vessel disease 
under stress. Areas of low CDE amplitude 
were found to match with the scar regions, 
whereas high CDE amplitude was found to 
correlate with areas of viable, ischemic 
cardiac tissue [Pesola, 1999], Kandori et al. 
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Table 1: Survey on several MCG studies in the field of CAD (for abbreviations see text). 


recorded 64-channel MCG in 8 patients 
with angina and 4 healthy subjects and 
analysed the relationship between currents 
recorded before and after exercise test as 
well as after interventional therapy. They 
calculated the exercise induced current for 
each channel as ratio of current after 
exercise to that during rest. CAD patients 
displayed three distinct patterns in current- 
ratio maps depending on the affected 
coronary vessel. After successful therapy 
the pattern returned to that of healthy 
subjects [Kandori, 2001], 

As a further source parameter 
equivalent current dipoles (ECD) were 
calculated on the basis of an iterative non¬ 
linear optimization algorithm at selected 
times during the QT interval. Dipole 
orientation in the xy plane as well as 
source strength were determined as the 
average of the calculated values over 5 ms 

for each time point and subject. While source strength did not differ between groups, dipole orientation deviated in most of the CAD patients from 
the normal values at at least one of the time points, in particular at T maximum and T offset (sensitivity 80 %, specificity 90 %) [van Leeuwen, 1997, 
1999], 





healthy 
n = 15 

Sin = 2.1 ±0.5 


61 channel system 

Fig. 1: Spatial distribution of QT interval over the thorax. 


CAD 


n =16 

Sin = 2.6±1.0 


PARAMETERS OF THE MAGNETIC FIELD 


Due to the detection and localization of abnormal current flows associated with ischemia, MCG reveals independent, complementary informa¬ 
tion compared to body surface ECG [Lant, 1990]. In order to register and quantify this additional information, the magnetic field was analyzed in 
course of the cardiac cycle using various parameters. As the current flow in the normal myocardium is almost uniform during ventricular repo¬ 
larization, the orientation is less uniform when ischemia is present. In one approach magnetic field orientation, defined as the angle between the line 
joining the centers of gravity of the positive and negative field components and the right-left line of the torso, was calculated during the QT interval 
at rest. These values separated CAD patients with normal ECG from the control group (sensitivity of 85%, specificity of 68% [van Leeuwen, 1999], 
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rest after bicycle exercise 

Fig. 2: Stress induced changes in orientation 
during ST segment (group mean). 


Orientation, quantified using the maximum spatial gradient of the magnetic field, was also used 
to determine changes in the ST segment and T wave during exercise (figure 2). In 27 patients with 
single vessel disease, Hiinninen et al. could show alterations in the magnetic field distribution of 
the mean ST-segment and T-wave apex during exercise testing [Hiinninen, 2000], Changes in the 
ST-segment orientation were most profound immediately after cessation of stress, whereas those in 
the T wave merged later. In a further study they assessed the optimal recording locations in 
magnetocardiographic mapping under stress in patients with specified ischemic regions referred to 
coronary angiography and no previous myocardial scar. Using magnetic isofield maps of the ST 
segment and the T wave amplitudes they could identify optimal locations for the detection of 
changes caused by ischemic injury current dependent on the affected coronary vessel. The 
identification of these regions could be optimized by calculating a discriminant index (DI) for each 
subgroup in each recording location [Hanninen, 2001], A further improvement of this approach was 
given by taken the heart rate into account. In 24 patients with single vessel disease Takala et al. 
could show that heart rate adjusted magnetic field map rotation over the ST segment and 20 ms 

interval on the T wave apex after exercise testing improved the detection of ischemia caused by stenosis > 75 % [Takala 2002], 

Beside the fact that the current flow becomes less uniformly oriented in the presence of ischemia, the associated magnetic field gets weaker. This 
results in a decrease in the maximum integral value in the iso-integral maps obtained for QRS and ST-T time intervals. Such changes have been 
reported by Tsukada et al. who showed a reduced ratio of maximum re- and depolarization values in ST-T and QRS iso-integral maxima compared to 

those of the control group [Tsukada, 2000], Sato et al. confirmed these data on 25 patients with stenoses > 75 %. 88 % of the patients showed either a 

rightward shift of current arrows, a multidipole pattern or decreased maximum time-integral values [Sato, 2001], 

Examining the information content of 
the magnetic field maps in a group of 15 
patients with suspected CAD, the 
Karhunen-Loeve Transformation (KLT) 
was used to calculate the first 12 
eigenvectors of each QRST integral map 
which in turn permitted the assessment of 
the non-dipolar content of the rest and 
stress integral maps. Under stress, the non- 
dipolar content revealed differences 
between patients with and without relevant 
stenoses [Van Leeuwen, 1999b], 



Fig. 3: Left: Example for CDV map category 0 and 4. Right: Relative number of CDV maps found 
in each of the categories before and after interventional therapy using boxplots; c= category. 


UNSHIELDED SYSTEMS 

Major limitations for widespread 

clinical use of MCG are the cost of the 
equipment including the need of 

electromagnetic shielding. With the development of devices not requiring shielding the availability of these systems in clinical settings should 
increase thus allowing the performance of MCG measurements in acute disease states. Fenici et al. could show repolarization abnormalities due to 
myocardial ischemia on the basis of magnetic field parameters first described by Hanninen et al [Fenici, 2002] [Hanninen, 2002], On the basis of 
current density reconstruction within the ST-T interval at rest differences could be shown between normals and a group of 52 CAD patients using a 
one channel unshielded system. These data based on a classification system derived from the quantification of the symmetries or asymmetries of the 
maps that will be influenced by disturbances in repolarization [Hailer, 2003], The results could be confirmed in a further study investigating patients 
in the course of coronary interventional therapy with a four channel magnetometer. While the CDV maps differed between CAD patients and healthy 
subjects before intervention, map distribution returned toward that of normals after successful therapy (figure 3) [Hailer, 2002], At present 
measurements were performed in more than 400 patients with suspected CAD. Consistent with the prior data CAD patients could be identified on the 
basis of CDV map classification with a sensitivity of 73,3 % and a specificity of 70,1 % [Hailer, 2004], 


CONCLUSION 


The future value of MCG in the detection of CAD will depend on the availability of these systems in a clinical environment, e.g. in chest pain 
units. The fact that MCG is not easily accessible for patients with acute coronary syndrome could explain why clinical data in this field are less 
extensive that in the area of arrhythmias. One may also question whether angiography as gold standard is appropriate in assessing patients with CAD. 
Newer work shows that patients with signs and symptoms of the disease may well demonstrate abnormal values in different clinical tests, including 
MCG [Van Leeuwen, this volume]. There is still a demand for fast and noninvasive screening methods for patients with suspected CAD or the follow 
up of patients after interventional therapy. A first step towards this goal was made by the development of systems with no need for shielding. Yet, the 
poorer quality of the signal compared to those obtained in a shielded location remains a disadvantage. A further improvement of this technique would 
include standardized data acquisition as well as data analysis procedures. This will result in a fast and generally accepted interpretation of the 
recordings, comparable to 12-lead ECG. The potential of MCG should be used for closing the still existing gap in the noninvasive diagnosis of CAD. 
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Magnetocardiogram is a useful noninvasive tool to identify patients at high risk of atrial fibrillation and 

heart failure hospitalization 
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Atrial fibrillation (AF), the most common sustained arrhythmia, is commonly due to reentrant mechanism. Conduction slowing in the atria 
predisposes to atrial reentry and AF. We hypothesized that magnetocardio-gram (MCG) which detects small difference in electrophysiological 
phenomena may be a sensitive tool to identify atrial conduction slowing and patients (pts) at risk of AF. 1) We compared P wave duration during 
sinus rhythm as determined by MCG and electrocardiogram (ECG) in 36 consecutive pts (18 with and 18 without prior history of AF). MCG P wave 
was significantly longer in pts with prior history of AF (99 ± 9 ms) than those without AF (80 ± 12 ms, P < 0.01) but there was no difference in ECG 
P wave between the pts with and without AF. 2) 22 pts with ischemic heart disease and no prior history of AF were divided into 2 groups: group 1, P 
wave on MCG < 115 ms (11 pts, 107 ± 5 ms) and group 2, P wave on MCG > 115 ms (llpts, 131 ± 13 ms) and were prospectively followed. There 
was no significant difference in ECG parameters between the 2 groups. Prolonged P wave duration by MCG identified those with low ejection 
fraction and ventricular dilatation. During the follow up (2.1 ± 0.7 years), group 2 pts had a higher incidence of AF (18 %) and heart failure (HF) 
hospitalization (45 %). Group 1 pts had no AF and a lower incidence of HF hospitalization (9 %). In conclusion, MCG is a useful tool in identifying 
pts with abnonnal cardiac hemodynamics and risk of adverse cardiac events, which can not be identified by ECG. 
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To innocuous electric median nerve stimulation at wrist, N20-P30 responses have been generally accepted as the earliest cortical responses (ECR). 
N17 activity is essentially subcortical in origin as it survives in patients with cortical or thalamic lesions abolishing N20-P30 and later cortical 
responses. Pre- and intra-operative recordings have shown a negligible contribution of N17 to cortical potentials. By recording magnetic fields 
generated by cerebral currents, magnetoencephalography (MEG) offers a non-invasive approach to evaluate sensory cortical functions in humans. 
Compared with EEG recordings, MEG reflect intracellular rather than volume currents, and magnetic signals are essentially transparent to extra¬ 
cerebral structures. Thus, MEG gives more accurate source localization and quantitative estimation of neuronal activation in cortical fissures than 
EEG. Using MEG recordings, the ECR at 20-30 ms (N20m-P30m) to electric median nerve stimulation, corresponding to N20-P30 responses by 
EEG, is generally localized in the area 3b of the primary somatosensory (SI) cortex. Earlier excision studies have shown that removal of SI cortex 
leads to the abolishment of the tactile ECR. However, the exact functional roles of the cortical generators for tactile ECR during perception of 
stimulus inputs remain not clear. 

In contrast to innocuous stimulation, noxious stimulation of the skin elicits cortical activation in bilateral second somatosensory (SII) cortices. As 
reported in earlier scalp and intra-cortical EEG recordings, the ECR peaking at 140-180 ms to C02 laser stimulation of the hand is generated in the 
cortex of the upper bank of the sylvian fissure, although SI cortex activation has been reported in some studies. MEG’s preferential sensitivity to 
fissural signals offers promising approach for non-invasive studies of human cortical pain processing. Although the SII cortex may play some roles in 
tactile learning and sensorimotor integration, its functional role in noxious infonnation processing is not clear. 

In the present talk, in addition to a concise introduction to the electromagnetic evidence and features of ECR elicited by innocuous electric as well as 
noxious stimulations, I will highlight some functional characteristics of underlying earliest cortical generators in processing non-painful and painful 
stimulations in humans. 


Finding the fast pathways through the visual system: MEG’s contributions 

J.M. Stephen 

Department of Radiology, University of New Mexico, Albuquerque, NM USA 

Monkey studies have shown that multiple parallel pathways carry infonnation through the visual system depending on visual characteristics and 
location of the stimuli in the visual field. The extension of this work to humans has been slow due to the difficulty of obtaining high spatial and 
temporal resolution, noninvasively. Although fMRI has helped map out many of the >20 visual areas previously identified in monkey studies, the 
limited temporal resolution has not allowed fMRI to fully describe the visual system. With more widespread access to sophisticated modeling 
approaches for MEG, the timing of different visual areas is providing significant results. Using circular sinusoids presented in central and peripheral 
locations (2.3 and 24 degrees, respectively), we found that dorsal stream visual areas did not have significantly different onset times, whereas ventral 
stream areas had progressively longer onset times as the visual information traveled from VI, V2/V3 to V4 [1], In addition, we found that 
peripherally placed stimuli had faster onset times in the intraparietal sulcus and superior lateral occipital gyrus (putative area MT) compared to 
centrally placed stimuli. These results agree with previous studies that showed equivalent onset times in VI and parietal occipital sulcus [2,3]. 
Ffytche et al. [4] found that motion stimuli with a fast drift rate activated V5 before activating area VI. These results suggest that multiple visual 
pathways exist in humans, that allow rapid transmission of critical visual infonnation to relevant cortical areas. In general, the fast transmission 
pathways connect dorsal visual stream areas. 

[1] Stephen, J.M., Aine, C.J., et al. 2002. Vision Research 42, 3059-3074. 

[2] Aine, C.J., Supek, S., George, J.S. 1995. Int J Neurosci 80, 79-104. 

[3] Vanni, S. Tanskanen, T„ Seppa, M„ Uutela, K., Hari, R. 2001. Proc Natl Acad Sci USA 98, 2776-2780. 

[4] ffytche, D.H., Guy, C.N., Zeki, S. 1995. Brain 118, 1375-1394. 
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Middle latency auditory evoked fields reflect dispersion effects in the cochlea 

A. Rupp 1 , N. Sieroka 1 ’ 2 , T. Dau 3 

'Section of Biomagnetism, Department of Neurology, Heidelberg University, 2 D-GESS, ETH Zurich, Switzerland, 3 Centre for Applied 

Hearing Research, Technical University of Denmark, Lyngby, Denmark 

Source analyses of middle latency auditory evoked fields (MAEF), using an optimal chirp stimulus (100-10,000 Hz) that compensates for the travel¬ 
time differences along the basilar membrane, show that phase delays between channels in the auditory pathway are preserved up to primary auditory 
cortex [1], Further, tone pulses (1100 Hz) embedded in complex tone maskers - linear raising (m+) and falling (m-) frequency sweeps from 200-2,000 
Hz and harmonics in zero-phase sine phase (mOsine) - produce different psychoacoustic detection thresholds, due to the dispersive characteristics of 
the cochlea [2], We investigated these dispersion effects and the correspondence between neuromagnetic responses, psychoacoustical loundness, and 
(simulated) auditory-nerve (AN) activity. The magnitude of the MAEF source-waveforms, when plotted against the temporal position of the tone 
pulse relative to the masker period, showed: (i) an U-shaped pattern for the m+ masker, (ii) the same pattern but smaller values for the m- masker, 
and (iii) an inversed U-shaped pattern for the zero-phase sine masker. The loudness scales and the simulated AN-activity revealed similar patterns. 
The close correspondence between the neuromagnetic, psychoacoustical and AN-activity patterns suggest that MAEFs might serve as an objective 
tool to investigate peripheral dispersion effects non-invasively. Thus, such data provide constraints on models of peripheral filtering in the human 
auditory system. 

[1] Rupp, A. et al. 2002. Hear. Res., 174, 19-31. 

[2] Kohlrausch, A., Sanders, A. 1995. J. Acoust. Soc. Am., 97, 1817-1829. 


Intracranial responses in auditory cortex of the monkey and human: 

Potential impact on interpretation of MEG 

Mitchell Steinschneider 1 and Matthew A. Howard 2 . 

Albert Einstein College of Medicine 1 and University of Iowa Medical Center 2 , USA 

Proper interpretation of neuromagnetic responses requires a detailed analysis of the underlying neural generators that can only be accomplished 
through intracranial recordings. Our studies have combined comprehensive intracranial recordings in the auditory cortex of the awake monkey with 
intracranial recordings in patients undergoing surgical evaluation for medically intractable epilepsy. Investigations in monkeys utilize multicontact 
electrodes that simultaneously record across laminae auditory evoked potentials (AEPs), the derived current source density (CSD), multiunit activity, 
and frequency-specific EEG changes. Studies in the human examine AEPs and EEG changes using intracranial depth and subdural grid electrodes. 
Results in both humans and monkeys indicate multiple overlapping source generators of neural activity in close proximity to each other that extend 
across auditory cortical fields. These neural generators, both within and across cortical fields, evoke different patterns of activity to the same stimuli. 
Higher frequency oscillations in the gamma range are often better indices of underlying neural organization than the lower frequency components 
present in the AEP Eaminar analysis of evoked activity in the monkey reveal complex and overlapping laminar patterns of both excitatory and 
inhibitory activity beginning in the thalamorecipient zone and extending over time to include supra- and infragranular laminae. We conclude that 
great caution must be exercised when attempting to extrapolate detailed features of auditory cortical organization from non-invasive recording 
techniques. 

Supported by DC00657 and DC042890. 
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Static and dynamic representation of complex sounds: from tonotopy to musical notes 

Takako Fujioka 1 ' 3 , Laurel J. Trainor, Bernhard Ross 4 , Ryusuke Kakigi' and Christo Pantev 4 
'The Rotman Research Institute, Baycrest Centre for Geriatric Care, University of Toronto, Canada 2 Departmant of Psychology, 
McMaster University, Canada, department of Integrative Physiology, National Institute for Physiological Sciences, Japan, institute 

for Biomagnetism and Biosignalanalysis, University of Munster, Germany 

When listening to speech or music, those complex sounds are encoded in the auditory system in two distinct ways of static and dynamic 
representations before being accessed by higher cognitive processes. Static infonnation of acoustic features such as the frequency spectrum can be 
indexed by auditory evoked responses (e.g. source strength and location of Nlm), whereas dynamic encoding can be observed in context-based 
responses (e.g. mismatch negativity (MMN)). The presentation will focus on our recent studies of the two ways of representation for complex sound 
in human auditory cortex. First, tonotopic organization of the Nlm generator is demonstrated up to 14kHz, and in a different manner for the complex 
tones for the superior-inferior direction. A more pronounced tonotopic representation (more distant location between center frequencies) in the right 
hemisphere gave evidence for right hemispheric dominance of spectral processing. Second, dynamic aspects of encoding complex sound were 
investigated using MMNm, the magnetic counterpart of the mismatch negativity. The infrequent melodic deviations in a single melody for pitch 
contour and interval structure, and even alterations in two-simultaneous melodies elicited MMNm suggesting that melodic information is 
automatically processed in auditory memory trace as an abstract code regardless the absolute pitch information. The significantly larger MMNm 
amplitudes in musicians suggests that musical training enhances the ability to automatically register changes in melodic structure. 


Sounds of one hand clapping 

Tommi Raij 

MGH/MIT/HMS Athinoula A. Martinos Center for Biomedical Imaging, Charlestown, MA, USA 

There are other ways than sounds to activate the auditory system. For example, omission of an expected sound will activate several auditory areas. 
Such "omission responses" have been shown in a wide range of species including humans. Random omissions (7%) of tone pips in an otherwise 
regular stimulus sequence evoked strong MEG responses over temporal and frontal areas, with peak amplitudes at 145-195 ms. Response amplitudes 
were 60% weaker when the subjects were not attending to the stimuli. Omission responses originated in supratemporal auditory cortices bilaterally, 
indicating that auditory cortex plays an important role in the brain's modelling of temporal characteristics of the auditory environment. Omission 
responses could be generated locally due to disruption of steady-state temporal activation patterns, and/or globally through long-range top-down 
effects from higher-order attention-dependent (e.g., frontal) areas. 

Also visual stimuli that have pre-leamed auditory associations (e.g., letters of the alphabet) can activate auditory areas. In an MEG study with 
auditory and visual letters, the auditory and visual brain activations first converged around 225 ms after stimulus onset and then interacted in the right 
temporo-occipito-parietal junction (280-345 ms) and the left (380-540 ms) and right (450-535 ms) superior temporal sulci. For these stimuli, the 
unimodal visual activations can also spread to supratemporal auditory cortices, probably through learning-dependent top-down connections from 
multisensory association cortices. 
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Specificity for pitch and other sound features in subdivisions of human auditory cortex? 

A. Gutschalk 

Eaton-Peabody Laboratory, Massachusetts Eye and Ear Infirmary, Boston, MA 02114 & Research Laboratory of Electronics, 

Massachusetts Institute of Technology, Cambridge, MA 02139, USA 

Specificity within the auditory cortex is not yet well established, and has only recently been shown for tonal pitch. This paper summarizes a series of 
MEG studies on tonal pitch using click trains. Periodic click trains are perceived with a well defined pitch at repetition rates above ~30 Hz. No pitch 
is perceived when clicks are repeated at the same average rate, but the spacing is jittered. Contrasting periodic and jittered click trains allows for the 
recording of auditory evoked fields (AEF) that are evoked by the pitch of the stimulus (the ‘pitch response’, PR). The PR consists of a subcomponent 
of the sustained field and the Nlm. Source analysis tracks both to the lateral part of Heschl’s gyrus, consistent with earlier PET and fMRI data by 
Griffiths, Patterson et al. The PR shows selective adaptation within the sequence; while the Nlm mainly depends on the interval between the 
repetition of related events, the sustained field adapts to the relative duration between related and unrelated events on a time scale of tens of seconds. 
The PR can be separated from the AEF that is equally evoked by periodic and jittered click trains in both its source location and specificity. The 
location of the non pitch sensitive AEF is more posterior, towards the planum temporale. The posterior sustained field varies strongly with the 
stimulus intensity. Moreover, the posterior Nlm is strongly attenuated when changes in the sound occur in the absence of intensity transients. 


Plastic change of auditory coding: 

Effect of hearing loss and developmental change 

Nagamine, T. ', Fujiki, N. 2 ’ 3 , Morita, T. 2 ’ 4 and Takeshita, K. L 5 

1 Human Brain Research Center, 2 Department of Otolaryngology-Head and Neck Surgery, Kyoto University Graduate School of 
Medicine, Kyoto, Japan; 3 Department of Otolaryngology and Broncho-Esophagology, Otsu Red Cross Hospital, Japan; 4 Brain 
Research Unit, Low Temperature Laboratory, Helsinki University of Technology, Finland; 5 Department of Pediatric Neurology, 

Omiya Medical Association Hospital, Japan 

Functional organization of auditory cortex has variable alterations depending on insults and developmental change. These changes can be manifested 
by the derangement of auditory evoked magnetic fields (AEFs) to various kinds of auditory stimuli. This paper summaries our studies of patients 
with deafness and healthy children. 

In patients with unilateral deafness, difference of N 100m latency between the hemispheres disappeared. No change just after the removal of 
peripheral part indicated this disappearance was produced in time range of months. As for the dipole strength, it was not affected by disease duration, 
except for cases of congenital deafness or early onset deafness who produced enhanced ipsilateral response. This might show possibility of 
disinhibition of ipsilateral pathway. 

Loudness sensation such as loudness recruitment was demonstrated by steeper gain increase in patients with inner-ear disturbance. Increment of 
dipole moment of NIOOm was more rapid according to the stimulus intensity in patients than that in healthy subjects. This was also shown 
transiently in patients with sudden deafness only when patients had subjective disturbance. 

In school-aged children, later component N250m could catch up with rapid rate stimulation such as 2-s interval, with which NIOOm could rarely be 
evoked. N250m was located more medial, anterior and inferior to NIOOm. In longer stimulus intervals, N250m was dampened and NIOOm was 
enhanced, showing functional discrepancy. 

Plastic change of auditory system can be shown by NIOOm in adults, although the temporal profile differs by its cause. However, N250m can be 
better serving as a functional benchmark in children. 
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Resolving multiple simultaneous sources of visual evoked magnetic fields 

Seppo P. Ahlfors 

MGH/HMS/MIT Athinoula A. Martinos Center for Biomedical Imaging, Charlestown, MA, USA 

Most visual stimuli are expected to evoked temporally overlapping activity in several functionally different cortical areas located close to each other. 
This is a special challenge to attempts of modeling the sources of visually evoked magnetic fields and scalp potentials. Combining information from 
different techniques, for example MEG and EEG, can help to resolve simultaneous sources that may be ambiguous in either one alone. The 
information provided by functional MRI can be used to guide the MEG/EEG source analysis by suggesting likely sources among the ones that can 
explain the recorded MEG/EEG data. Subspace regularization [1] provides one theoretical framework for incorporating fMRI information. Analysis 
of the sensitivity patterns of MEG/EEG sensors as well as cancellation of signals due to multiple simultaneous sources is of importance when 
combining multimodal data. The interpretation of results should take into account that some activation may be missed because of limitations of the 
recordings techniques. One straightforward by potentially very useful way how fMRI can help interpret MEG source estimates is by identifying the 
side of a sulcus a or a gyrus where the source is likely to reside. This can also resolve whether the source current is orientated inwards or outwards 
with respect to the cortical surface. On the other hand, the orientation of the source current can typically be detected accurately; a change in the 
orientation can suggests the presence of multiple neural populations, even when they may look like a single extended region in fMRI. Examples of 
multimodal recordings in individual subjects demonstrate the usefulness of combing multiple types of data to better determine the patterns of cortical 
activity related to visual information processing. Supported by the Whitaker Foundation and the MIND Institute. 

[1] Ahlfors, S.R, Simpson, G.V. 2004. Geometrical interpretation of fMRI-guided MEG/EEG inverse estimates. Neuroimage 22, 323-332. 


MEG and ERP correlates of visual motion and object-related attention 

Mircea Ariel Schoenfeld 

Dept, of Neurology II, University of Magdeburg, Germany 

Studies have shown that visual attention modulates the processing of motion. Event-related potentials (ERP) and magnetic fields (ERF) were 
combined functional magnetic resonance tomography (fMRI) to precisely analyze these modulations with regard to their temporal characteristics and 
corresponding sources. Attending coherently moving dots elicited an effect in the timerange 100-350 ms which was maximal over occipito-temporal 
electrode/sensor sites. The electric/magnetic fields of this effect were modeled with two sources each located next to the left and right temporo- 
occipital junction. Their location was coincident with the hemodynamic activations in the fMRI. These findings were used in a second study to 
investigate the timing and cortical localization of feature integration mechanisms in object-based attention. Subjects attended to one of two 
superimposed transparent surfaces formed by arrays of dots moving in opposite directions. The spatio-temporal analysis revealed evidence for a rapid 
increase in neural activity localized to a color-selective region of the fusiform gyrus when the surface moving in the attended direction displayed an 
irrelevant color feature. The results provide support for the “integrated competition” model of object-selective attention and point to a dynamic 
neural substrate for the rapid binding process that links relevant and irrelevant features to form a unified perceptual object. 
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Neural mechanisms of attentional focusing in visual search 

Jens-Max Hopf 

Department of Neurology II, University of Magdeburg, Germany 

Event-related electric potentials (ERP) and magnetic fields (ERMF) were employed to investigate the neural mechanisms that underlie the attentional 
resolution of ambiguous feature coding in visual search. We compared neural activity related to target discrimination under conditions of major 
versus minor feature overlap between target and distractor items and found that increasing feature overlap lead to a focal enhancement of neural 
activity in ventral occipito-temporal areas, consistent with the larger need to attenuate distractor interference. These results suggest that distractor 
attenuation proceeds as a stepwise operation, with different spatial locations containing interfering features being suppressed successively, and 
support theories of visual search that emphasize location-based attentional selection as a key mechanism in resolving ambiguous feature coding. 
Further we investigated whether attentional selection of nonspatial features operates in a location-independent manner. A visual search task was used 
in which the spatial distribution of nontarget items with attended feature values was varied independently of the location of the target. The presence 
of task-relevant features in a given location led to a change in ERP/ERMF activity beginning ~140 msec after stimulus onset, with a neural origin in 
the ventral occipito-temporal cortex, that was independent of the location of the actual target. This effect was followed by lateralized activity 
reflecting the allocation of attention to the location of the target, which began at ~170 msec poststimulus. Source analyses indicated that the 
allocation of attention to the location of the target originated in more anterior regions of occipito-temporal cortex than the feature-related effects. 
These findings suggest that target detection in visual search begins with the detection of task-relevant features, then allowing spatial attention to be 
allocated to the location of a likely target, which in turn allows the target to be positively identified. 


Temporal Dynamics of Attention/Working Memory 

C.J. Aine 1 ’ 2 

'Albuquerque VA Medical Center, 2 University of New Mexico SOM, Albuquerque, NM USA 
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Magnetic Imaging of Action Currents and Potentials in reveal the bidomain properties in cardiac tissue 

F. Baudenbacher, J.R. Holzer, V. Y. Sidorov, L. Fong 
Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA 

The origin of the magnetocardiogram (MCG) and the relative information content between the MCG and the electrocardiogram (ECG) remain central 
questions in biomagnetism. To provide key insights to this question, we have recorded spatially resolved MagnetoCardioGrams (MCGs) and 
optically imaged the transmembrane potentials generated by planar wave fronts on the surface of the left ventricular wall of Langendorf perfused 
isolated rabbit hearts. The MCGs were combined and analyzed to produce a time series of two-dimensional action current maps. Overlaying epi- 
flourescent images allowed us to identify a net current along the wave front and perpendicular to gradients in the transmembrane potential. This is in 
contrast to a traditional double layer model where the net current flows along the gradient in the transmembrane potential. Our findings are supported 
by numerical simulations treating cardiac tissue as a bidomain with unequal anisotropies in the intra - and extracellular spaces. Our measurements 
reveal for the first time the anisotropic bidomain nature of cardiac tissue during plane wave propagation. These bidomain effects play an important 
role in the generation of the MagnetoCardioGram and potentially lead to new information. 
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SQUID vs. Optically Pumped Magnetometer: a comparison of system performance 

G. Bison 1 ,A. Pasquarelli 2 , A. Weis' and S.N. Erne 2 

'Physics Department, University of Fribourg, Switzerland 
2 Biosignals and Imaging Technologies Department, ZIBMT, University of Ulm, Germany 


ABSTRACT 

Due to their extremely weak intensity, biomagnetic fields are usually detected by means of superconducting interference devices (SQUID) which 
are currently the most sensitive magnetic sensors available. To reduce construction and running costs, extensive efforts have been dedicated to the 
development of High-Tc SQUIDs working at liquid nitrogen temperature and suitable for biomagnetic applications, but the results up to now are only 
partially satisfactory. 

Optically pumped magnetometer (OPM) technology, although known for many years, has recently emerged as a promising alternative to 
SQUIDs, at least for applications with more moderate signal-to-noise requirements. Since OPMs were presented to the biomedical community, a 
controversy has surrounded their proposed usefulness for different biomagnetic applications. This contribution compares SQUID and OPM methods 
based on their physical principles: for example by the fact that a SQUID measure the flux through the pick-up-coil area, while the OPM is sensitive 
to a volume averaged magnetic field component. At the implementation level, the most significant differences between SQUID and OPM are, 1) the 
absence in OPMs of cryogenic and vacuum parts reducing the gap between source and sensor, and 2) the need for more complex read-out electronics 
of the OPM. 

The intrinsic insensitivity against RF-interference of the OPM reduces the shielding requirements and may ease the introduction of biomagnetic 
systems into noisy environments. Full technological maturity is by far not yet achieved, thus leaving a large improvement potential, which may 
eventually allow the OPM to compete with the SQUID even in high-resolution applications. 

KEY WORDS 

Instrumentation, SQUID-magnetometer, optically pumped magnetometer, biomagnetic field detection. 

INTRODUCTION 

SQUID magnetometers are to date the most sensitive devices available for the detection of magnetic fields. Their extraordinary sensitivity allows 
to pick-up extremely weak magnetic signals like the ones associated to the bioelectric currents in living organisms. These signals have different 
intensities ranging in orders of magnitude from hundreds of pico-Tesla for magnetocardiogram (MCG) down to units or even fractions of femto-Tesla 
for magnetoneurogram (MNG) and brain-stem associated activities. 

To exploit at best the SQUID’s sensitivity, several techniques were adopted to deal with the external sources of magnetic and electromagnetic 
disturbances, which surpass the biomagnetic signals by orders of magnitude. Typically available techniques are passive and/or active shielding, 
gradiometric design, filters, data-processing algorithms and so on, all having pros and contras especially under the point of view of costs/benefits 
ratio. Yet SQUIDs are superconducting devices, which need an appropriate cryogenic technique, they are usually operated in liquid helium, and this 
requires maintenance at regular intervals and at relevant costs. Furthermore a SQUID system cannot be switched-on and operated in minutes like 
every ECG recorder off-the-shelf, here the start-up time from cool-down to full operation takes several hours. Therefore two ways of operations are 
usual: 1) if sporadic measurements have to be carried-out, they need to be carefully scheduled and possibly grouped to take the most benefit from 
cooling-down the sensor for a limited period of time, and then the system warns up, 2) for routine measurements the system is kept cold and running 
continuously for several months. 

All these arguments may have strongly limited over the years the wide spreading of the biomagnetic method. To overcome these problems many 
efforts have been done in searching for less expensive technologies, both for the sensors and for the noise suppression techniques. The discovery of 
high-Tc superconductivity in the late 80’s looked very promising since it allows one to operate sensors with liquid nitrogen at much lower costs and 
system complexity, but to date the reliability of HTc-SQUIDs is still an open issue. 

Recently it was shown that OPMs can be used for the detection of MCG. The reported noise level was more than one order of magnitude larger 
than that of currently available SQUIDs, but this is just the beginning of the study and there are good margins for improvements. This paper presents 
a comparison SQUID vs. OPM and an estimation of chances for the future. 


METHODS 

A DC-SQUID consists of a superconducting ring interrupted by two identical Josephson junctions. If a biasing current is applied to the ring, 
under normal conditions the current will equally split into the two junctions and the phase of the wave functions will be the same. An external 
magnetic flux threading the ring determines a phase shift of the two wave functions, which combined with the flux quantization phenomenon leads to 
the well known periodic flux-to-voltage transfer function, whose period is equal to a flux quantum 4 >q of 2 ■ 10‘ 7 Gauss cnr. Beside the extra¬ 
ordinary sensitivity, this behavior has the disadvantages of nonlinearity and limited dynamic range. Therefore the SQUID is driven by suitable 
feedback-electronics that keeps the working point of the SQUID locked at a fixed value (Flux-Locked-Loop) [Clarke, 1969], [Koch, 1989], Such a 
SQUID magnetometer shows an excellent linearity over five to six orders of magnitude and the minimum detectable flux can be in the range of |I(|>q 

A typical biomagnetic SQUID-based system consists of: sensor array, cryogenic probe and readout electronics at room temperature. The 
cryogenic probe is required because the SQUID operates in liquid Helium at the temperature of 4.2 K, in a thermally shielded enclosure (dewar). The 
outstanding features of these instruments include a wide dynamic range (at least 100 dB), a wide bandwidth from DC up to a few tens of kHz at a 
noise level in the white region (Nyquist), ranging from a guaranteed values below 10 fT/VHz in commercial systems, down to 2 fT/VHz for the best 
research devices [Drang, 1995], 

In an OPM a vapor of paramagnetic atoms sealed in a glass cell (Fig. 1 b) is used to detect magnetic fields via the induced Larmor precession of a 
long-lived spin polarization [Alexandrov, 1993], Its principle of operation is based on a resonant circularly polarized laser beam, which creates spin 
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polarization in the atomic vapour via optical pumping. The magnetization precesses around the external magnetic field B 0 at the Larmor frequency 
W L which is proportional to the modulus of B 0 . A reference phase is imposed by a magnetic field B\ (radio-frequency) co-rotating with the spins 
around B 0 at frequency U) r{ . The driven magnetization induces an oscillating component of magnetization on the laser-beam axis, which modifies the 
optical absorption properties of the medium, which is detected by monitoring, with a photodiode and a lock-in amplifier synchronized to W rf , the 
beam power transmitted through the sample. The signal shows a resonant behavior when 6U lf approaches the Larmor frequency. The sensitivity of the 
Mx magnetometer to external magnetic fields varies as sin 26, where 0 is the angle between the laser-beam direction and the magnetic field and has 
thus an optimum sensitivity for 0 = 45. Therefore the bias field B 0 = 5 |JT is best applied to the sensor at 0 = 45° (45° geometry, see Fig. l.b). The 
photocurrent is analyzed by a lock-in amplifier referenced by the radio frequency U) l{ . When operated in this Mx configuration the device acts as a 
scalar magnetometer whose signal is a function of |5 0 |. As the MCG-induced field changes B are much smaller than the offset field B 0 applied in the 
z direction, the signal changes of the magnetometer will be (to first order) proportional to Bz only; this makes the device act as an effective vector 
magnetometer. Active noise suppression, based on a first-order gradiometric design, can be arranged with two identical magnetometers. 



Figure 1 a) Geometry of a SQUID system, b) 45° geometry of the optical Mx magnetometer, c) Flat geometry of future optical magnetometers. 


The typical intrinsic sensitivity achieved with a cylindrical glass cell (diameter 20 mm, length 20mm), filled with a droplet of caesium metal and 
a buffer gas mixture, is around 100 fT/VHz [Bison, 2003], Recently the system could be optimized to operate with an intrinsic sensitivity of 63 
fr/VHz and a detection bandwidth of 140 Hz. Simple model calculations comparing the signals measured by SQUID and by OPM assuming a 
dipolar source shows that the absence of a dewar, i.e., a smaller distance between source and sensor, can partially compensate for the smaller 
sensitivity of the OPM. The effect is of course small for deep lying sources (an improvement of the sensitivity of only few percent) but can be 
relevant for sources near the surface (up to 40% for a depth of 30 mm). Future designs of optical magnetometers, based on a slightly different 
magnetooptical technique, will increase this effect significantly by using flat geometriey (Fig. 1 c). 


DISCUSSION 

The experience gained until now from the recording of MCG data in a weakly shielded or unshielded environment is showing that the intrinsic 
noise limits of the measurement systems (dominated in SQUID systems by the environmental noise) do not differ by orders of magnitude, but only 
by a factor 3 or 4. There is founded expectation that in the near future that gap can be reduced perhaps to a factor of 2. The advantage of the OPM 
based system, being free of cryogenics and practically using standard electronic equipment is clearly in favor of a broad clinical application of the 
OPM technique. The most important step will be to demonstrate that in the clinical routine there will be accepted applications where the slightly 
reduced sensitivity is not a decisive drawback. 
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Magnetoencephalography with an atomic magnetometer* 

M.V. Romalis, A. Ben-Amar Baranga, D.E. Hoffman 
Physics Department, Princeton University, USA 

Recently developed spin-exchange relaxation free (SERF) atomic magnetometers have magnetic field sensitivity competitive with low-temperature 
SQUID magnetometers. They also have high spatial resolution and can be easily expanded to multiple channels [1]. We are presently constructing a 
256-channel atomic magnetometer system suitable for human magnetoencephalography. The magnetometer consists of a glass cell containing a vapor 
of K atoms and two lasers used for optical pumping and probing of spin precession. Using a 2-D photodiode array the magnetic fields will be 
simultaneously measured on a two-dimensional grid in a cubical cell approximately 7 cm on the side located about 2 cm away from a human head. 
For repetitive signals it will also be possible to construct a 3-D magnetic map by scanning the pump laser, increasing the effective number of 
channels to several thousand. The high density of measurement points will allow detailed studies of localized brain activities, such as evoked by 
somatosensory or auditory stimuli. Numerical simulations indicate that spatial localization of the magnetic field sources within the brain should be 
improved by more than an order of magnitude compared with traditional SQUID systems. Since the magnetometer does not need a cryogenic dewar, 
it can be placed in a compact human-size magnetic shield with a high shielding factor and made entirely from non-conductive materials to reduce 
thermal magnetic noise. We have completed construction of most magnetometer components and will report on the technical details, numerical 
simulations and latest test results. 

Supported by NIH and Packard Foundation 

[1] I. K. Kominis, T. W. Komack, J. C. Allred and M. V. Romalis, Nature, 422, 596 (2003). 
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Argos 500: Operation of a Helmet Vector-MEG 
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ABSTRACT 

After four years of work a magnetoencephalographic system has become operational at the University of Ulm. The sensor design consists of a 
163 vector-magnetometers array oriented and located in a suitable way to cover the whole head of the patient. Four additional triplets are available as 
references to build software gradiometers. The helmet shaped sensor system is positioned to accommodate a lying patient in a fashion similar to MRI 
and CT scanners. Simultaneously to the MEG, 32 EEG channels (expandable to up to 128) and other relevant patient’s information up to a total 
number of 660 acquisition channels can be recorded. Sampling rate is fixed at 8200 Flz. All basic pre-processing, band-pass filtering, decimation and 
noise compensation are perfonned on the fly using an array of DSP. 

The vector-magnetometer array can be used as is or alternatively to emulate a classical array of magnetometers (synthetic magnetometers), with 
the choice of a variable geometry. In any measurement site it is possible, on the basis of data from the vector-magnetometer triplet, to calculate the 
component of the signal normal to the surface of the volume conductor, reducing the effect of the volume currents. 

The sensor is placed in a magnetically shielded room consisting of an external 8mm aluminum layer for eddy-current shielding and three layers of 
1.57 mm thick soft magnetic material. 

The global system is optimized for clinical use: the system requires only one step for maintenance per week. Preliminary measurements on 
normal subjects are presented. 
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INTRODUCTION 

The MEG system installed at the University of Ulm is specifically designed for clinical application and routine use, to allow investigation of a 
large number of patients per day. To reach this goal, the system design meets the requirements of reliability, high field sensitivity, minimal set-up 
overhead before each measurement and an easy-to-handle user interface. 

METHODS 

The MEG sensor system consists of a large vertical dewar, held on a fixed 
basement, with a side helmet-shaped cavity to accommodate the patient’s head 
during the measurement. All sensors (integrated SQUID magnetometers) are 
arranged in triplets to build vector-magnetometers. SQUID’s specifications include a 
square pick-up coil (12 mm in diagonal), integrated Applied Positive Feedback 
(APF) circuit and feedback-coil, intrinsic noise spectral density ranging between 2 
and 5 fT/VHz in the white region. 492 basic sensors (166 vector-magnetometers) 
cover the helmet-shaped sensing area and detect the magnetic fields of interest. 12 
additional sensors (4 vector-magnetometers), conveniently placed at 4 different 
locations around the helmet, are dedicated to the detection of the noise field. These 
two sensors arrays can be software combined to form synthetic vector-gradiometers. 

Moreover the vector-magnetometer array can be alternatively used to emulate a 
classical array of magnetometers (synthetic magnetometers), with the choice of a 
variable geometry. In any measurement site it is possible, on the basis of data from 
the vector-magnetometer triplet and of anatomical data from the patient, to calculate 
the component of the signal normal to the surface of the volume conductor, reducing 
the effect of the volume currents. 

The patient handling set-up is made of non magnetic material and responds to 
the same operating criteria of other imaging techniques like CT or MR scanners with 
the patient lying on the bed and being slide into the measuring apparatus. 

The sensor system is operated inside a magnetically shielded room (MSR) 

[Pasquarelli, 2000], with one aluminium layer eddy current shielding for AC stray 
fields and three layers of soft high-permeability shielding material for the lower 
frequency range. The inner room dimensions are 2.9, 2.9, 2.3 m (w, 1, h) and require 
a placement of dewar and patient’s bed along the diagonal of the room to ease the 
access for medical personnel and patients in the case of measurements and also for 
the technical staff and equipment in the case of maintenance. An overview is given 
in Fig. 1. 

All SQUID-drive-boards are mounted on top of the cryostat, arranged in 11 groups of 8 boards, each one driving 6 SQUIDs. All operations are 
remotely programmed by software, through fiber-optics link and dedicated decoder, communicating with a proprietary noiseless protocol. 

Up to additional 128 channels are available for EEG other physiological signals like respiration, temperature etc. and for triggers and parameters 
from other devices like stimulators, timers, supervisor devices etc. 

All conditioned signals are fed, trough adequate EMI filters, to the A/DC unit outside the MSR. The data conversion runs at 8200 Hz per channel 
and has 20 bits of resolution. The digital outputs of over 600 A/D Converters are collected, formatted in 32-bits words and sent to the acquisition 
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Fig. 1 - Overview of Argos 500 inside the MSR. 
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computer by a single chip device working in a star-like architecture. At this point the data stream is fed over an optical gigabit network link to the on¬ 
line DSP processing unit, for compensation (synthetic gradiometers), filtering and decimation to a final sampling rate which can be defined between 
62 and 4100 Hz according to the measurements requirements. The last steps in the chain are storage on the data server and archiving on removable 
media. 

In the perspective of clinical routine use of the system, significant efforts were made for the implementation of the acquisition and the post¬ 
processing systems, in order to reduce the global measurement time and the analysis time for patient. A side product of the DSP work is a slower data 
stream for the on-line monitoring of the signals. Here two different continuous displays are available: one low-resolution monitor with several small 
traces, which can show all signals or just selected groups among several possible choices, and one high-resolution monitor which shows a selection 
of up to 16 traces with richness of details. This second workstation is at the same time the operator’s interface for management and control of the 
measurements. The processing console provides a flexible and fast tool for data analysis, with both time domain and magnetic source imaging 
methods. 


RESULTS 

First of all the system was calibrated to obtain characterization data for both 
mechanical parameters and sensitivity of all sensors (this work is presented 
elsewhere in this conference) [Pasquarelli, 2004], 

After this process also the compensation matrix was calculated and the 
synthetic-gradiometers could be built. At this point it was possible to measure the 
background noise of the system with and without compensation. Sample spectra 
are shown in Fig. 2. 

The ultimate check was performed with real measurements of neurological 
signals, and specifically SEF elicited by median nerve stimulation at the wrist, 
and voluntary motor activity if the index-finger. In both cases clear signals where 
detected (see Fig. 3 and Fig 4). 



Fig. 2. Noise spectra of one sample channel without 
(black) and with compensation (gray). The white 
noise level shows just minor changes, whereas the 
1/f region is improved by 1 order of magnitude. 




Fig 3. Somatosensory evoked field elicited from median 
nerve stimulation at the wrist. Two traces with opposite 
polarities are shown. 


Fig. 4. Motor activity detected after voluntary movements of the 
index-finger of the right hand (101 Avg., 3 minutes recording). Two 
traces with opposite polarity are shown. 
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ABSTRACT 

Development of the CTF MEG system has been advanced with the introduction of a computer processing cluster between the data acquisition 
electronics and the host computer. The advent of fast processors, memory, and network interfaces has made this innovation feasible for large data 
streams at high sampling rates. We have implemented tasks including anti-alias filter, sample rate decimation, higher gradient balancing, crosstalk 
correction, and optional filters with a cluster consisting of 4 dual Intel Xeon processors operating on up to 275 channel MEG systems at 12 kHz 
sample rate. The architecture is expandable with additional processors to implement advanced processing tasks which may include e.g. continuous 
head localization/motion correction, optional display filters, coherence calculations, or real time synthetic channels (via beamfomier). 

We also describe an electronics configuration upgrade to provide operator console access to the peripheral interface features such as analog signal 
and trigger I/O. This allows remote location of the acoustically noisy electronics cabinet and fitting of the cabinet with doors for improved EMI 
shielding. Finally, we present the latest performance results available for the CTF 275 channel MEG system including an unshielded SEF (median 
nerve electrical stimulation) measurement enhanced by application of an adaptive beamformer technique (SAM) which allows recognition of the 
nominal 20 ms response in the unaveraged signal. 
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INTRODUCTION 

There is increasing demand for real time processing of MEG data. The CTF MEG electronics was originally designed only for 3 real time 
processing tasks: anti-alias filter, sample rate decimation, and higher order gradiometer formation. However with the added requirements for 
continuous head localization, crosstalk correction, and advanced tasks such as real time beamformer output, the hardware architecture is no longer 
able to meet the processing demands. The bottleneck is not in the individual processing boards but in the backplane interconnects between boards. As 
an alternative to design of a new hardware architecture, it was determined that a computer processing cluster could meet the present demands and 

offer better flexibility for future requirements. A computer cluster is a network of fast PCs connected together with a high speed network. Many of 

the world’s fastest supercomputers are now based on these clusters. Utilization of clusters for processing tasks has only recently become feasible 
because the important components are becoming sufficiently fast: processors (Intel 3 GHz Xeon processor - 3Gflops/s), memory (18Gbytes/s cache 
and 3Gbytes/s RAM), and networks (Myrinet - 180 Mbytes/s). The implementation of the cluster in the CTF MEG system is described in the 
following section along with additional improvements to the user’s peripheral equipment interface. 

Preliminary performance data for the CTF 275 channel MEG system has been reported previously [Fife, 2002], Additional noise measurements 
are reported with more details in the ‘performance’ section below. In the same section, we present an additional example of the unshielded system 
performance via a somatosensory evoked field (SEF) measurement. 

ELECTRONICS 

The CTF MEG electronics architecture has been described in some detail previously [Fife, 2002], The basic architecture is unchanged in the 
present advancement, however a computer processing cluster has been interposed between the control bank (consisting of DSP/PGA processing units 

and a data collection control unit) and the acquisition computer. This is shown 
in the block diagram of fig. 1 where the inputs from SQUID sensor and EEG 
electrode amplifiers to the MEG and EEG banks have been omitted to simplify 
the diagram. Output from the stimulus computer to the various stimulation 
devices (e.g. auditory, visual, and somatosensory) has also been omitted in the 
figure. The cluster processing block resides within the electronics cabinet while 
the user’s peripheral interface (UPI) has been moved from the cabinet to the 
operator’s console where the acquisition and stimulus computers are located. 

The cluster, as initially configured, consists of 4 nodes (off the shelf rack 
servers) and is expandable as additional processing tasks are identified. One 
computer acts as the server node (N 0 ) and houses the fiber optic interface to the 
electronics control bank while the other 3 are slave nodes (Ni - N 3 ). Each node 
has dual Xeon 3 GHz processors running Redhat 9, 512 Mb of level 2 cache, 
1024 Mb of RAM, and 40 Gb hard drive. The nodes are interconnected using a 
high speed network interface (Myrinet, Myricom Inc., Arcadia, CA, USA) for 
data and gigabit ethemet in parallel for housekeeping. The acquisition computer 
and each node has a Myrinet interface card connected to an 8 port Myrinet switch unit via fiber-optic cable. 

Tasks presently implemented in the cluster, in order of process flow, are: a primary input buffer (accepts data from the control bank), crosstalk 
correction, filter/decimation, gradient balancing, optional filtering, and output buffer (sends data to the acquisition computer). The input buffer and 
part of the crosstalk task run on N 0 , remaining crosstalk task and filter/decimation run on Ni, gradient balancing and optional filtering run on N 2 , and 
the output buffer runs on N 3 . 

The central challenge for the implementation of the cluster at 12 kHz is the process timing. When processing one sample at a time, each task must 
execute the input-process-output cycle in 83.3 ps. However, it is known that the Myrinet I/O processing has a point of maximum efficiency (64KByte 
transfers) and that task processing incurs a time penalty due to a "cache hit" where the task program and data must be loaded from RAM into the 
CPU cache. Tests have shown that processing multiple samples at the same time spreads the cache hit penalty and the EO overhead across each 



Figure 1. Electronics block diagram 
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sample in the group to lower the overall sample execution time when 
compared to single sample processing. For example, timing tests run on the 
gradient balancing engine revealed a single sample processing time of about 
24 ps while the simultaneous processing of 11 samples resulted in about 19 
ps per sample. 

The remote UPI provides convenient access to trigger and analog I/O 
connections as well as remote master power/reset buttons. It also allows the 
noisy electronics cabinet to be located at some distance from the operator. 

Features of the redesigned UPI include 16 ADCs, 4 DACs, a video time 
channel, system time channel (for synchronization), and 16 general purpose 
TTL digital channels configurable as either input or output. The DAC 
channels can be programmed to generate a variety of waveforms. The 
analog data is filtered and decimated in exactly the same manner as MEG or 
EEG channels. 

PERFORMANCE 

The CTF 275 channel MEG system perfonnance has been reported 
previously [Fife, 2002] as a preliminary noise spectral plot that compared 
shielded and unshielded response. Additional data collected in a typical 
shielded room demonstrated average white noise of about 4 fT rms/VHz 
while low frequency noise remained below 10 fT rms/VHz for frequencies 
higher than about 0.2 Hz. An examination of individual channel noise 
spectra revealed that, in the frequency range from 1 to 2 Hz, all but 3 
channels had noise between 3.5 and 7 fT rms/VHz (data collected at 
National Institute of Mental Health, Bethesda, MD). 

An additional illustration of system perfonnance is the response to median nerve SEF in fig. 2. This measurement was done in an unshielded 
environment with 628 trials averaged from data collected at a sample rate of 1200 Hz with DC to 300 Hz bandwidth. Note the stable baseline in the 

System performance improvement by increased channel count was 
demonstrated by enhanced beamformer source resolution [Vrba, 2004], We 
use a minimum variance beamformer (SAM) which provides a spatial filter 
based on a point dipole forward model [Robinson, 1999], Virtual sensors 
(VS) are formed by applying beamformer weights for a selected 3D 
location in the brain to the MEG measurement. An example of SAM 
performance comparison with standard MEG channel output using the 275 
channel system is provided in fig. 3 for the SEF data of fig. 2. The left hand 
side of the figure compares a segment of unaveraged data over 4 trials from 
the MEG channel with largest SEF signal with a VS calculated at a location 
of maximum SEF response. Comparison of the data with the trigger marks 
reveals clearly resolved single trial SEF response from the VS which are 
not resolved in the MEG sensor. A single trial is isolated in the right hand 
side of the figure for a comparison of the unaveraged VS response with the 
averaged VS and averaged MEG sensor response. The single trial VS 
clearly resolves the expected general SEF morphology while the averaged 
response compares favorably with the averaged MEG sensor response 
morphology, but with better signal to noise ratio. 

CONCLUSIONS 

We have described the introduction of a computer processing cluster in the electronics of the CTF MEG system. The present configuration with 4 
nodes is capable of anti-alias filtering, sample rate decimation, higher gradient formation, crosstalk correction, and real time head localization at 
sample rates up to 12 kHz for a total of 464 channels (up to 304 MEQ 128 EEQ and 32 stimulus). Future implementations with additional nodes 
should make it possible to include advanced real time processing tasks such as SAM virtual sensors and user supplied algorithms. 

The perfonnance of the CTF 275 channel MEG system has been demonstrated with an example of a median nerve SEF paradigm. The system 
provides good DC response and allows the detection of the 20 msec component from an unshielded single trial via SAM virtual sensor. 
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Figure 3. SEF signal comparison of MEG channel output with 
SAM virtual sensor output. Vertical lines are stimulus trigger 
points. 
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Figure 2. Averaged SEF data example; top: projection plot of 
individual channel data, bottom: all channels overlaid. 
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Active Shielding Method for an MEG device 

J. Simola, S. Taulu, L. Parkkonen, and M Kajola 
Elekta Neuromag Oy, Helsinki, Finland 

A method for external interference suppression based on a novel dual feedback arrangement is described. An MEG device with the individual 
channels run in the flux-locked-loop is provided with a second negative feedback loop using the sensors of the array itself as zero indicators, and a set 
of coils around the array as actuators. The set of coils, typically three to six, is tailored to optimally compensate for the dominating external 
interference field patterns, and linear combinations of the sensors, optimally sampling the interference at the very location of the array, are chosen for 
zero indicators. The number of control quantities (linear combinations) can be smaller or equal to the number of compensation coils. No separate set 
of dedicated reference sensors is required. Both the external interference and the compensation provided by the actuator coils are signals from 
sources outside of the sensor array, so that their effect can be separated from the biomagnetic signals by the Signal Space Separation (SSS) method 
[1], This leads to recovery of the biomagnetic signal with no distortion even in those channels used in the zero indicators. By this method an 
interference suppression by a factor exceeding 50 has been achieved. In this way even magnetometers can be run with a relatively light magnetic 
shield because the achieved interference suppression is similar to that of reasonably balanced gradiometers. Owing to this the magnetometers have a 
dynamic range similar to that of gradiometers, as far as external interference is concerned. However, contrary to gradiometers, the biomagnetic 
signals, even from deep sources, are not attenuated. This compensation method can be used both with and without a magnetically shielding room. 

[1] Taulu S., Kajola M., Simola J. 2003. The Signal Space Separation method. Proc. of the 14 th Conference of ISBET (Nov. 2003, Santa Fe, NM) 
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Signals in Micro tesla Magnetic Fields 
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ABSTRACT 

We have built and demonstrated a system with open geometry that measures biomagnetic signals such as magnetoencephalogram (MEG), 
magnetocardiogram (MCG) and magnetomyogram (MMG) simultaneously with ultra-low field nuclear magnetic resonance (NMR) free induction 
decay signals (FID). The system employs SQUID gradiometers and operates in measurement fields of 2-50 DT, with proton Lannor frequency is in 
80-2000Hz range. A pre-polarizing field in the range of 2-30 mT is generated by a resistive coil. Two different types of SQUID gradiometers were 
used: an axial second-order gradiometer and a tangential thin-film planar first-order gradiometer. The gradiometers were placed inside a fiberglass 
dewar at about 1 cm distance from a subject. All measurements were performed inside a single-layer magnetic shielded room. 

KEY WORDS: Low field NMR, MRI, FID, MEG, MCG, MMG, SQUID, gradiometer, Biomagnetism. 

INTRODUCTION 

In order to understand biomagnetic data, it must often be combined with anatomical data obtained by MRI. SQUID sensors used to acquire 
biomagnetic signals can not operate in large magnetic fields required for conventional MRI, consequently two completely separate systems must be 
used to acquire the biomagnetic and anatomical data. The possibility of recording NMR signals in small magnetic field opens up the possibility of 
acquiring tomographic images simultaneously with biomagnetic signals, as the magnetic fields required for imaging are now compatible with SQUID 
sensors [Volegov, 2004], NMR line widths are expected to be significantly narrower at very low fields. For a fixed relative inhomogeneity, 
broadening of the NMR line scales linearly with the strength of the measurement field providing NMR line-widths that approch the lifetime limit and 
raise the possibility of MRI at lower gradients [McDermott, 2002], Susceptibility artifacts causing spurious field gradients in the sample that broaden 
resonance lines are significantly reduced at low fields [Macovsky, 1993], These effects combine to offer the possibility that imaging at low fields may 
provide a regime of high sensitivity and resolution. Low field MRI may ultimately lead to direct imaging of the magnetic fields from bioelectric (for 
instance, neural) currents by observation of change of FID or echo parameters caused by local biomagnetic fields [Xiong, 2003], For such studies, 
simultaneous biomagnetic measurements will allow us to compare observed MR responses with a direct electrophysiological measure. It may be 
useful to employ MRI for current localization while using MEG to precisely track bioelectrical currents. 

This paper describes a system with open geometry that we have built and demonstrated to measure biomagnetic signals such as 
magnetoencephalogram (MEG) [Volegov, 2004], magnetocardiogram (MCG) and magnetomyogram (MMG) [Espy, 2004] simultaneously with ultra- 
low field NMR signals. We have obtained 1H NMR spectra from a variety of samples including water, mineral oil, and live tissue, and gradient 
encoded data from a water plastic phantom was used to construct simple 2D projection images [Matlachov, 2004], 

INSTRUMENTATION 

The system consists of a SQUID-based gradiometer, a prepolarizing coil and a coil system to generate uniform measurement field and gradients. 
Two different gradiometer configurations were tested: a second-order axial wire-wound gradiometer and a first-order tangential planar gradiometer. 
The gradiometer was placed inside a fiberglass liquid helium dewar model BMD-5 from BTi. 

The second order gradiometer consists of one bottom, one top and two center 14 mm diameter turns with 50 mm baseline. Two center turns are 
wound in opposite direction with respect to both bottom and top turns to create the second order gradiometer configuration with a balance level of 
about 0.1%. It has 4.4nT/<I>o field sensitivity and about 10fT/VHz resolution. The tangential gradiometer is Nb thin-film device deposited on one 
20x60 mm 2 silicon substrate composed of two 20x20mm" pick-up coils with 40mm baseline and balance level better than 0.01% [Stolz, 1999], The 
field sensitivity is 0.28nT/<J> 0 and field resolution better than 8fT/VHz. The pre-polarizing field, B P , was generated by a variety of different wire- 
wound coils of varying size and shape that were used in different experiments. For example, two 65mm diameter coils with 250 turns each were used 
for many phantom experiments. The two coils were positioned co-axial to the tail of the cryostat one coil right on the end of the tail and second coil 
about 60mm above. The coupling of B p with the axial gradiometer was decreased by more than a factor 10 4 by adjusting the relative coil positions. 
Only one pre-polarizing coil can be used with the tangential gradiometer, however the geometry is such that very little field is coupled. Rotating the 
gradiometer around the coil axis can further minimize mutual coupling. For simultaneous biomagnetic and NMR signals detection, an 80mm ID, 
160mm OD pre-polarizing coil, with ~800 turns was used. This coil could generate 30mT fields. A square Helmholtz-pair coil with 56cm sides was 
used to generate the measurement field, B m , orthogonal to B P . The tangential gradiometer was oriented orthogonal to both B m and B P coils to 
minimize flux coupled into the SQUID. We also used a 1.2m diameter round Helmholtz coil pair for simultaneous MCG and NMR detection The tail 
of the dewar containing the gradiometer was placed directly over the subject, minimizing standoff. The whole system is placed inside a single layer 
magnetically shielded room consisting of an aluminium eddy current shell and mu-metal static field shield [Sullivan, 1989], The shielding factor of 
the MSR is about 40 at 1 Hz and about 200 at 60 Hz, and the residual DC magnetic field is below 0.1 pT. 

METHODS 

The typical measurement protocol consisted of the pre-polarization of the subject tissue for a time that is comparable with the expected Ti (about 
2 s for water and about 1 s for living tissue). While B P is applied the SQUID is effectively turned off. The measurement field B m is always on. At time 
t = 0, B P was slewed to zero in less than 0.5ms. The SQUID electronics are activated with a time delay from 2 ms to 20 ms after B P reaches zero and 
signal acquisition was begun . We typically used about 1 ms delay after SQUID electronics turned on to trigger the data acquisition system, although 
the recovery time of the electronics was only 20ps. 

The electronics time delay could not be set below 2ms because of presence of large transient signal that suppressed the SQUID transfer function. 
The origin of this transient needs further investigation. This transient signal lasts up to 100ms with a few nT amplitude. A double exponential 
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approximation of the transient was used to digitally reduce the artifact in the data by more than a factor of 1000. Accurate approximation of this 
effect is very important for simultaneous MEG and FID signal detection, although it is less important for only FID detection. We speculate that this 
problem can also be eliminated by using an echo technique. 

We tested both axial and tangential gradiometers for simultaneous biomagnetic and NMR measurements. An axial gradiometer has the con¬ 
venience that B m can be oriented in any direction in perpendicular to gradiometer axis plane, however axial gradiometers detect the maximum 
biomagnetic signal when it is placed offset with respect to the primary bioelectric source. A tangential gradiometer detects the maximum 
biomagnetic signal when it placed immediately above the source so it receives both biomagnetic and FID signals from the same area. A tangential 
gradiometer, however, must be oriented parallel to both B m and bioelectric current source, which may be inconvenient. Using two orthogonal 
tangential gradiometers placed at the same point, one can orient the first gradiometer to detect only biomagnetic signal and the second gradiometer to 
detect only FID signals. In this case B m should be oriented perpendicular to bioelectric source. 


RESULTS 

We have experimentally verified, using a water phantom and a 65mm ID pre¬ 
polarizing coil, that 80% of the NMR signal originates from the top 20mm of the 
sample, while 20% of the signal strength originates from more distant regions. This 
leads us to conclude that it is unlikely to observe a change of FID or echo parameters 
caused by bioelectric activity because we receive only a small portion of signal from 
the area with strong primary bioelectric currents. 

A two-coil pre-polarization technique was designed to improve sensitivity for 
deeper sources. We used two different sizes and polarity of coils to minimize the FID 
signal originating from the proximity of the gradiometer. One 150mm ID coil 
(17G/A) was placed in the plane of the dewar tail surface while a smaller 65mm ID 
coil (42G/A), was placed 38mm above the first one. Figure 1 shows the experimental 
data for FID amplitude vs. distance between the dewar tail and the surface of a 
400ml beaker of water. Zero distance corresponds to water touching the tail. The 
three sensitivity curves correspond to pre polarization with the large coil only, small 
coil only and both coils. Figure 1 clearly illustrates that when only a single coil is 
used, the majority of the FID signal arises from the upper water layer. When two 
coils with opposite polarity were used the top-most ~15mm of water contributes 
minimally to the FID signal (FID signal doesn’t change vs. distance), and most of the 
FID signal arises from deeper layers. This fact was also verified using a phantom 
with two fluid layers having different T 2 times. 

High signal to noise ratio (SNR) of the FID signal is an important requirement 
for creating high quality MRI. Figure 2 shows the experimentally measured SNR for 
100 epochs recorded in real time from a water phantom and from human subject at 
B. = 20 mT and B m =10 pT. The real time SNR is around 5 for water and around 2 
for living tissue. Two deep minima at 48-th and 84-th epochs for the water curve due 
to trapped flux jumps in the SQUID. SNR can be improved by averaging as well as 
by increasing B P . Averaging 100 epochs increases the FID signal SNR to about 20 
for human subjects. SNR can be improved by another factor of 10 by increasing B P 
to 200mT (which can be realized in practice). 

Both long term and short term stability of the instrument are also important. Our 
instrument showed about 0.1 Hz frequency deviations at 390Hz during one 10- 
minute run and 0.15 Hz deviation between two runs separated by 30 minutes. 
Stability of other FID parameters such as T 2 time and phase jitter needs further 
experimental investigation. 
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of Energy, Office of Biological and Environmental Research. 
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Successful Neurorehabilitation Therapy May Recruit Dysfunctional Brain Areas 

T. Elbert, D. Bulach, G. Urswatte, M. Meinzer, E. Taub, B. Rockstroh 

University of Konstanz, Germany and University of Alabama, Birmingham, USA 

Cerebrovascular accident (CVA) is the leading cause of disability in Europe and in the United States. The consequent motor deficits that patients 
sustain result in significant reduction of the quality of life and also in very large costs. The reduction of CVA-related disability thus represents a 
health care priority. 

Edward Taub and coworkers [1] have developed a set of techniques involving intensive, extended massed-practice that randomized controlled studies 
indicate can substantially reduce the motor deficit of patients with mild to moderately severe chronic strokes. For the arm, the therapy involves 
promoting use of the more-affected upper extremity by: 1) most importantly, intensive training of that extremity with the behavioral technique termed 
shaping for several hours per day over a period of consecutive weeks (massed practice), and 2) constraining use of the less-affected ann for a target 
of 90% of waking hours with a padded mitt that prevents compensatory use of the hand. We have used conventional magnetic source imaging of 
representational cortex and in addition a new MEG-based method for mapping functional and dysfunctional brain areas after stroke (ASWAM) to 
further elucidate the plastic brain processes associated with the administration of effective neurorehabilitation therapy. Results show that successful 
neurorehabilitation produces a large activity-dependent plastic reorganization in the brain that doubles the size of the portion of primary motor cortex 
involved in generating movements of the more-affected extremity. The mapping of abnormal magnetic brain waves indicates that dysfunctional 
neural networks have become modified in response to the treatment. 

Research was supported by the Deutsche Forschungsgemeinschaft. 

[1] Taub, E., Urswatte, G., Elbert, T. 2002. New treatments in neurorehabilitation founded on basic research. Nature Reviews Neuroscience 3, 228- 
236. 
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Brain sensorimotor hand area organization and functionality in acute stroke: insights from MEG 
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Universita Campus Biomedico, Rome, Italy. 


ABSTRACT 

A ‘snapshot’ of the acute phase following a cerebral stroke is fundamental to understand the brain’s adaptation potential, which plays a key role in 
clinical recovery. In the present study, the characteristics of the spontaneous cerebral activity of the rolandic region and the evoked response 
following the stimulation of the median nerve have been studied trough magnetoencephalographic (MEG) recordings in 32 patients admitted to 
neurological ward for first-ever acute ischaemic stroke involving the upper limb. 

Both in affected (AH) and un-affected hemispheres (UH) spontaneous activity evaluation was performed by analysis of spectral properties, 
individual alpha frequency and spectral entropy. Evoked activity was evaluated by characterization of the cortical sources activated by stimulation of 
the median nerve, in terms both of absolute values in each hemisphere, and of interhemispheric differences. Interhemispheric waveshapes cross¬ 
correlations were also carried out. Neurophisiological findings were correlated with neuroradiological ones and clinical scores. 

In the acute phase after an ischemic attack, the rest activity showed signs suggesting: a reduction of spectral richness; an increase of the intra- 
regional neural synchrony; an increase of the lower with respect to the higher frequency powers. 

Signs of enhanced excitability were present in the AH following a cortical lesion, usually in combination with preserved hand functionality. An 
enhanced excitability of the UH was paired with larger lesions with cortical involvement. 

KEY WORDS: Sensorimotor cortex, hand representation, spontaneous activity, neural plasticity, interhemispheric asymmetries, middle 
cerebral artery, magnetoencephalography (MEG). 

INTRODUCTION 

Studies of short- and long-term changes in neural circuitry and brain connectivity after an acute ischemic brain lesion represent the basis of our 
still partial understanding of the post-stroke recovery. Functionality of sensorimotor brain hand areas were studied in the early post-stroke stages both 
for understanding acute neuronal reaction properties and posing the basis for studying the variation between acute and stabilized conditions. The 
“resting state” characteristics of the neuronal activity are of primary importance in determining the brain’s processing capabilities. In particular 
characterization of different functionality levels of specific areas have also been assessed in the rest state, as shown extensively for alpha activity 
asymmetries between the dominant and the non-dominant hemisphere [Cobb 1963]. On the other side, rest and activation properties are not directly 
related; this is mainly due to the fact that different neuronal populations contribute to the recorded signal in the two conditions. 

We have studied through magnetoencephalography (MEG) the organization of the neural activity in the cortical areas devoted to sensorimotor 
hand control at rest and under electrical stimulation in patients who had suffered an infarction in the middle cerebral artery (MCA) territory within 
the previous 10 days. 

METHODS 

Thirty-two patients (mean age 68 ± 12 years, range 30-86, 15 females, 17 males; 15 right and 17 left hemispheric stroke) were enrolled between 2 
and 10 days (mean 5.2 ± 2.6 days) after first-ever acute ischemic stroke affecting the arm-hand control diagnosed on the basis of clinical history and 
examination and confirmed by brain MRI. The criteria for inclusion were: clinical evidence of a motor and/or a sensory deficit of the upper limb and 
a neuroradiological diagnosis of ischemic brain damage. The local Ethical Committee approved the experimental protocol and written informed 
consent was obtained. Clinical state was assessed by National Institute of Health scale (NIH), Bartel index (BI) and 5 degree motor score (0-1.5, 0 no 
symptoms). 

MEG recordings were performed by a 28-channel system covering a total scalp area of about 180 cm' inside a magnetically shielded room 
(Vacuumschulze). Spontaneous activity from rolandic left and right regions was recorded (3 minutes-open eyes). 

Somatosensory evoked fields (SEF) were recorded following separate electrical stimulation of left and right median nerves at wrist via surface 
disks by 0.2 ms electric pulses (cathode proximal), with an interstimulus interval of 641 ms. Stimulus intensities were set just above threshold 
inducing a painless thumb twitch. Data were continuously acquired for 3 minutes (1000 Hz sampling rate, 0.48-250 Hz bandwidth) and processed 
off-line. 

Power Spectrum Density (PSD) was obtained at rest state via Bartlett procedure (0.5 Hz resolution, Hamming window, 60% overlap). Band 
power was calculated in delta (2-3.5 Hz), theta (4-7 Hz), alpha (7.5-12.5 Hz), betal (13-23 Hz), beta2 (23.5-33 Hz) and gamma (33.5-45 Hz) ranges; 
spectral entropy and the frequency of maximal value in alpha band (IAF) were also calculated. 

Equivalent current dipole (ECD) strength of the two initial waves (M20, M30) were compared with a nonnative data set obtained from a large 
population of sex-age matched healthy subjects [Tecchio 1997]. Waveshape analysis was carried out via a correlation coefficient comparing the SEF 
morphology in the two hemispheres [Tecchio 2000] on two different post stimulus time epochs: 30 ms starting from M20 onset, and the following 50 
ms. These two epochs roughly provide interhemispheric correlation of sensory processing taking place in primary sensorimotor areas (morph_Sl) and 
in parietal brain areas ( morph_par). 

RESULTS 

Rest activity. Cause no interhemispheric “physiological” differences in single bands were evidenced in controls, hemispheric comparisons were 
performed on values in AH with respect the UH. Results are summarized in figure 1. Low rhythms absolute power (delta and theta) were higher in 
AH with respect to the UH both for right and left side involvement. If delta absolute power in the two hemispheres were compared to control values, 
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control subjects 


not only AH were significantly increased (univariate ANOVA, [F(l, 39)=11.643, p= .002]) but also UH [F(l, 39)=4.609, p= .038], No effect was 
revealed in alpha and beta band. 

Only a reduction in gamma AH absolute powers were observed with respect to controls in right-damages patients (p=.012; p<.001, p=.002 
respectively). The symmetry found in the spectral entropy values in controls resulted lost in patients, with a decrease of entropy in AH (entropy in 
AH lower than in UH, p=.002). The same alteration was found for AH IAF, lower in AH than in UH (p< .001). 

Evoked activity. Eight patients out of 32 showed no identifiable SEF from their AH. These patients all presented severe clinical deficits, as 
testified by the poor sensory scores and NIH and BI values. Besides those cases with no response, six more showed a reduced absolute ECD strength 
in the AH. Two patients, both with parietal cortical lesions, showed an increased ECD strength in the AH; both had good levels of motor and sensory 
functionality. The remaining three patients showed an ECD strength increase in the UH with respect to the AH and belonged to a group of six patients 
with absolute ECD strength higher than normative values in the UH. There was a significant correlation between the dimension of ischemic lesion 
and the strength of the M20 component in the UH (rho = 0.43; p = 0.019). All patients with a UH response enhancement, had a cortical involvement. 

Rest/Evoked activity vs clinical picture. The following relationship with clinical scores were found: 

1. Correlation with AH gamma band power, showing higher 
gamma values with better clinical pictures (NIH Spearman’ 
rho = -.624, p=.001; BI rho = .637, p<001; motor score rho = 

.484, p=.010); 

2. Correlation with AH spectral entropy, indicating better clinical 
pictures paired with higher entropy values (NIH rho = -.430, 
p=.028; BI rho=.375, p=.059). 

3. Correlation with AH IAF values, indicating better clinical 
pictures paired with higher IAF values (NIH rho = -.440, 
p=.031; BI rho=.434, p=.034; motor score rho=.433, p=.035) 

4. Lower than normal absolute values of ECD strengths in AH 
resulted statistically correlated with more severe clinical 
pictures (NIH rho [32]=-.57; p = 0.001 and rho = -0.40; p = 

0.023 respectively for M20 and M30; BI rho [32]= 0.60; p < 

0.001 and rho [32] = 0.44; p = 0.013 respectively for M20 and 
M30; motor score rho = 0.52; p = 0.002 and rho = 0.37; p = 

0.037 respectively for M20 and M30; sensory score rho = 

0.35; p = 0.051 and rho = 0.30; p = 0.097 respectively for 
M20 and M30). 

5. Patients with abnormal inter-hemispheric asymmetry of the 
ECD M20 and M30 strength showed a worse clinical picture. 

6. Patients showing morphSI outside the nonnative range or no 
SEF response showed a significantly worse clinical picture. 



left-damaged patients 


right-damaged patients 



DISCUSSION 


Figure 1. Mean band powers in UH (fine line) and AH (bold line)for 
patients affected by left and right lesion. Asterisks indicate values 
statistically different in AH with respect to UH. 


A lesion affecting the MCA territory of one hemisphere, induces in the acute phase an increase of the AH rolandic areas power at rest and a 
reduction of spectral richness (lower entropy), suggesting higher synchrony of regional neuronal activity. The worst clinical pictures were associated 
to reduced AH gamma activity and to an unbalance of the neural networks generating the main rhythmic alpha activity with an increase of the lower 
rather than the higher frequency powers. The slight positive relationship between higher entropy levels and better clinical pictures suggests that 
lesions not affecting neuronal activity frequency richness allow the preservation of area functionality. A worst clinical picture was also related to the 
decrease of individual alpha frequency in the affected hemisphere, which had been also previously observed [Pfurtscheller 1986], [Giaquinto 1994], 
[Juhasz 1997], This indicates that an unbalancing of the neural networks generating the main rhythmic alpha activity, with an increase of lower with 
respect to the higher frequency powers, has a negative incidence on hand control functionality. 

Signs of enhanced excitability were present in the affected hemisphere following a cortical lesion, usually in combination with preserved hand 
functionality. An enhanced excitability of the unaffected hemisphere was paired with larger lesions with cortical involvement; signs compatible with 
an abnormal transcallosal transmission and intracortical function of inhibitory GABAergic inter-neurons in the AH were found subtending UH 
enhancement. Spared responsiveness from Brodmann’s area (BA) 2 and posterior parietal areas despite an altered response from BA 3b was found in 
six patients, combined to high hand functionality. Excessive interhemispheric differences of the investigated SEF parameters represent a potential 
probe for tracking brain deficiencies following a monohemispheric lesion. 

Present results in acute phase increase the knowledge of the mechanisms governing brain adaptation/reaction capabilities, for future efforts to 
establish therapeutic and rehabilitative procedures. 

REFERENCES 


Cobb WA, The normal adult EEG In: D. Hill and G. Parr, editors. Electroencephalography. New York: Macmillan; 1963. P. 232-249. 

Juhasz C, Kamondi A, Szirmai I. Spectral EEG analysis following hemispheric stroke: evidences of transhemispheric diaschisis. Acta Neurol Scand. 
1997; 96:397-400. 

Giaquinto S, Cohianchi A, Macera F, Nolfe G. EEG recordings in the course of recovery from stroke. Stroke. 1994; 25:2204-9. 

Pfurtscheller G. Rolandic mu rhythms and assessment of cerebral functions. Am J EEG Technol. 1986; 26:19-32. 

Tecchio F, Rossini PM., Pizzella V., Cassetta E., Romani GL. Spatial properties and interhemispheric differences of the sensory hand cortical 
representation: a neuromagnetic study. Brain Res. 1997; 767: 100-108. 

Tecchio F, Pasqualetti P, Pizzella V, Romani GL, and Rossini PM. Morphology of Somatosensory Evoked Fields: inter-hemispheric similarity as a 
parameter for physiological and pathological neural connectivity. Neurosci Lett. 2000; 287: 203-206 


43 







S10 


Spatiotemporal patterns of movement-related fields in stroke patients 
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ABSTRACT 

This study shows the analyses of brain functions between stroke patients and normal subjects by observing neuromagnetic fields during button 
pressing tasks. The measurement system included force measurement, visual stimuli presentation and MEG measurement. A 122-channel whole-head 
MEG system (Neuromag 122) was used in the study. A total of 18 subjects (11 post-stroke and 7 control subjects) participated in the study. The 
equivalent current dipoles for MRFs were estimated with the following parameters, i.e., goodness of fit (>75%), confidence volume (< 6000mm 3 ), 
and the duration of dipole (> 10msec). In addition to the single dipole analysis, Minimum Current Estimates was applied for source estimation since 
neural activities for stroke patients were observed at the motor cortex of the contralateral side as well as other areas of the brain. As a result, 
contralateral motor cortex was activated for the normal subjects at 50msec prior to the force onset, whereas ipsilateral motor cortex was activated for 
some stroke patients and patients’ dipole moments differed not only in their locations but also in their latencies, ranging from -150 to 50msec to the 
force onset. The results were in agreement with the findings by PET and f-MRI studies, therefore, it was implied that the compensating motor 
functions were shifted neighboring areas of the brain due to the recovering motor function after stroke. 

KEY WORDS 

Movement-related fields, minimum current estimates, motor cortical area, stroke, hemiparesis, reorganization, plasticity, motor recovery, motor 
control, dipole analysis. 

INTRODUCTION 

Stroke is a major disease to cause numbers of neurological dysfunctions including aphasia and paresis. Although many patients continually 
suffered from profound disabilities, some patients obtained variable functional improvement regardless of the initial severity. However, the 
mechanism of recovery is unclear so that there is no method to predict when and how the recovery is progressed for particular patient’s case. Several 
functional imaging techniques such as fMRI, MEG, PET, and TMS were tried to visualize activity in sensori-motor brain areas [Rossini, 1998] 
[Trompto, 2000] [Foltys, 2003], and to identify the mechanism of the reorganization of the damaged cortical areas. In particular, MEG demonstrates 
spatial localization for the active neuronal pools with an extremely high time resolution. Although there are still few approaches using MEG, recent 
study used MEG to show cortical area reorganization in stroke along with other functional imaging techniques [Rossini, 1998], They used current 
dipole method and showed the site of post-stroke region expressed by the location of current dipoles overlapped with MR images. They discussed 
that the concordant reshaping of temporally-sensitive functional imaging such as MEG and TCS should be used to interpret findings by fMRI. At the 
same time, they also indicated that the major drawback by using MEG is that MEG gives three-dimensional identification of sources obtained by 
means of an inverse procedure relying on the choice of a mathematical head model. In this study, minimum current estimates (MCEs), an estimate for 
selecting a good model structure for multi-dipole models [Uutela, 1999], were used to visualize MEG data into functional mapping for identifying 
reorganizing cortical areas after stroke by comparing with the images by normal subjects. 

METHODS 

A total of 18 subjects (11 post-stroke and 7 normal subjects) participated in the study. All normal subjects were male, right handed and their ages 
ranged from 22 to 25 years. All patients (ages 55-79, 7 males and 4 females, right handed) were affected by right or left hemiparesis. Six patients 
were also diagnosed as dysarthia. Average NIH stroke scales at the onset of stroke was 5.6 points. The subjects were instructed to press a pad with 
his/her right thumb continually at irregular intervals of three to 10 seconds with instructional visual prompts showed at the centre of the display. Two 
tasks were investigated; Simple pressing task and force-controlled pressing task. In simple pressing task, the intensity of applied force was not 
designated so that the subjects determined each force level individually. In force-controlled pressing task, additional visual infonnation indicating the 
intensity of force exerted was displayed as well as the range of the targetting force intensity, so that the subjects tried to exert their force within the 
intensity range. 

The measurement system included force measurement, visual stimuli presentation and MEG measurement. A 122-channel whole-head MEG 
system (Neuromag-122) was used in the study. The responses were recorded with a 0.03-100Hz band-pass filter (BPF) and digitized at a sampling 
rate of 400Hz. Averaging rejection was applied to the MEG data beyond 3000fT/cm. FSR measured the pressure of the thumb, and the top of FSR 
was covered with the rubber sheet for noise attenuation. MEG and FSR data were averaged off-line after the elimination of artifacts. The 
measurement data were averaged more than fifty times with the trigger signal of the intensity of the applied force of the thumb. The averaged data 
ranged from -1000 to 1000 milliseconds around the FSR onset time. The data were also processed with a BPF (0.1~30Hz) and SSP (Signal Space 
Projection) as noise rejecters. 

Based on averaged MEG signals, the amplitude of magnetic fields were determined for each scalp location by adding independently directed two 
differentiated spatial vectors given by a pair of orthogonally arranged SQUID sensors. The technique for determining the amplitude data for 
orthogonally arranged SQUID sensors can be found in [Sakaguchi, 2001], The amplitude data were categorized into five regions (left-frontal, right- 
frontal, left-motor, right-motor, and occipital). 

MEG data were further analyzed with two methods, the equivalent current dipole method and MCEs. Analysis using MCEs was detailed in 
[Uutela, 1999], In brief, MCE visualizes multiple sources by using Ll-nonn estimates, where explicit apriori information was not required. In this 
study, estimated source by MCEs was matched with the results of the equivalent current dipole method. The equivalent current dipoles for 
movement-related fields were estimated with the following parameters, i.e., goodness of fit (> 90%), confidence volume (<250mm 3 ), and the 
duration of a dipole (> 10msec). We estimated dipoles at 2.5msec intervals from the force onset (0msec). A spherical head model was used in the 
dipole fitting. Dipoles of the subjects were also matched with respective head anatomies using their MRI data to identify the source location of the 
signals geographically. The analysis was aimed to clarify the relationship between spatio-temporal changes in the intensity of the applied force by the 
thumb and the active intensity of the motor and somatosensory cortex. 
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Figure 1. MCEs obtained at between -150~50ms of the movement onset for normal 
subjects.(right hand movement) 


RESULTS 

Among patients group, eight patients 
had difficulty in pressing buttons enough 
to control their force levels, resulting in 
lacking the amplitude of trigger signals 
for averaging MEG signals in some 
experimental conditions. Three patients 
completed all trials. In order to complete 
two groups, three patients’ data were 
further analyzed. Figure 1 shows results 
of MCEs obtained by a normal subject. 

The tone-bar represents the intensity of 
source strength. After the movement 
onset, an increased activation in the 
contralateral motor area is apparent. This 
tendency was found in all normal 
subjects. The current dipole analysis 
confirmed that a dipole was located at the 
motor cortex on the contralateral hemi¬ 
sphere. Figure 2 is a typical example of 
MCEs obtained by a patient. In this case, 
affected hemisphere is on the left. An 
increased activation in multiple motor 
areas is apparent on the ipsilateral motor 
cortex, whereas relatively low activity can 
be seen on the contralateral hemisphere in 
spite of the time to the movement onset. 

DISCUSSION 

MCEs for patients showed that acti¬ 
vated areas were widespread compared to 
those for normal subjects before and after 
the movement onset. The spatial patterns 
of activated areas showed differently 
between patients, however, three locations 
were mainly activated; Anterior part of 

precentral gyrus on the unaffected hemisphere, posterior part of the precentral gyrus on the affected hemisphere, and supplementary motor area on 
the affected hemisphere. Activation on the unaffected motor cortex by post-stroke patients was observed by several fMRI studies [Cao, 
1998][Cramer, 1997], Such phenomena have been explained by the association of the anterior corticospinal tract. Our results implied that activation 
at the other two areas on the affected hemisphere may be due to the achievement of stroke recovery by supplying blood flow from neighboring areas 
that were not affected by the stroke. 
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Figure 2. MCEs obtained at between -150~50ms of the movement onset for the subject A 
with right hemiparesis.(right hand movement) 
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Intensive language training enhances brain plasticity in chronic aphasia 

MF Meinzer, TR Elbert, C Wienbruch, D Djundja, B Rockstroh 
Department of clinical psychology, University of Konstanz, PO-Box 25, 78457 Konstanz, Germany 

Focal clusters of slow wave activity in the delta frequency range (1-4 Hz) as measured by magnetencephalography (MEG) are usually located in the 
vicinity of structural damage of the brain. In the present study we investigated the change in Delta Dipole Density (DDD) in 28 patients with chronic 
aphasia (> 12 months after stroke) following cerebrovascular stroke of the left hemisphere before and after intensive speech and language training. 
Evaluation criteria included neurolinguistic testing (Aachen Aphasia Test Battery), brain imaging was performed using MEQ coregistered with 
structural magnetic resonance scans (MRI). 

The 28 patients received 30 hours of intensive language training over a period of 2 weeks (3 hrs./day). Neuropsychological testing and MEG 
assessment were performed before and after therapy. Patients improved on average after therapy and single case analysis revealed significant 
improvements in about 85% of the patients. Perilesional dysfunctional slow wave activity decreased in about 60% of patients that improved in the 
AAT, which might be an indicator of re-recruitment of these areas. An increase of perilesional activity in 9 patients with improved language functions 
might be explained by further segregation of dysfunctional components from the speech network with correspondent slow wave synchronisation. 
Changes in global measures of aphasia severity (AAT profile score & Token Test) varied with magnitude of change only in left hemispheric DDD. 
These results emphasize the impact of perilesional areas in the rehabilitation of aphasia even years after the stroke and might reflect reorganisation of 
the language network that provides the basis for improved language functions after intensive after intensive training. 


The influence of brain tumor treatment on pathological activity in MEG 

J.C. de Murick", A. de Jongli 1 ', J.C. Baayen 2 ’, R.M. Heethaar 1 W.P. Vandertop 2 ’, C.J. Stam 3) 

VU University medical center, Amsterdam, The Netherlands, Department of 
1) Clinical Physics, MEG center; 2) Neurosurgery; 3) Clinical Neurophysiology 

Interictal spikes represent only a very small part of the MEG/EEG data. It has been proposed in the literature to use another pronounced property of 
the MEG/EEG data of patients with epilepsy: the enhanced delta band activity. It was investigated whether pathological delta activity can be 
localized using simple dipole models, and in how much this delta activity is influenced by tumour treatment (neurosurgery and/or radiotherapy). 

MEG and MR recordings were made both before and after the treatment in 17 patients. Patients with malignant tumors exhibited more delta activity 
than patients with low graded tumors and meningiomas. In all patients with high grade tumors the delta power was lower after the treatment. 
Preoperatively 14 clusters of equivalent dipole sources describing focal activity were found in 12 out of 17 patients. Thirteen of these clusters were 
located near the tumor, and one cluster near an edema border. After the operation 13 such clusters were found in the same 12 patients (one exception). 
Eleven clusters were located near the lesion border and 1 cluster near the edema border. The shift in source locations after operation was in general 
considerably smaller than the dimension of the preoperative tumors. This finding indicates that similar areas generate the pre- and postoperative delta 
activity. Furthermore, focal delta sources were found in a case without tumor recurrence and also in cases that most tumor tissue was removed. These 
findings suggest that the pathology underlying the slow waves is the irreversible structural damage to the surrounding white/grey matter. 
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ABSTRACT 

Recently we developed a new method, intra-epidermal electrical stimulation (ES), for selective stimulation of the cutaneous A-delta fibers. In the 
present study, we recorded cortical responses to ES using magnetoencephalography to investigate the nociceptive processing in the human cortex. ES 
was applied to the dorsum of the left hand with a current intensity of 0.1-0.3 mA, which produced well-defined pricking sensations. The earliest 
cortical response was evoked in the primary somatosensory cortex (SI) at around 80-100 ms after the onset on the stimulus. Following the SI activity, 
the dorsal part of the anterior/mid insula and the secondary somatosensory cortex (SII) were activated almost simultaneously, whose activity were 
larger than that of SI in amplitude and peaked at around 150 ms. At a latency period of 200-300 ms, two deep sources in the anterior cingulate cortex 
(ACC) and anterior medial temporal area (MT) were active. The present results together with previous anatomical findings in animal studies suggest 
that SI and SII are serially activated following noxious stimulation but activations in SII and insula are parallel. Based on the activation timing in this 
study and projections from the insula to ACC and MT, we considered that ACC and MT are driven by insula. Therefore there are two nociceptive 
pathways: thalamus-SI-SII and thalamus-dorsal insula-ACC and MT. VPL and VMpo (Craig) are possible candidates of the thalamic nucleus 
responsible for the respective pathway. 

KEY WORDS 

MEG, insula, pain, SI, SII 

INTRODUCTION 

Although many recent studies using experimental animals and humans indicate that multiple cortical areas are involved in nociception, little is 
known about temporal profile of each cortical activity, especially in humans. Recently, we developed intra-epidermal electrical stimulation (ES) 
method [Inui, 2002a] for selective activation of cutaneous A-delta fibers. The peripheral conduction velocity of ES-evoked signals was about 15 m/s 
[Inui, 2002a,b]. This is an electrical stimulation, and therefore activates skin nociceptors directly without time delay unlike laser stimulation, which is 
suitable to investigate precise timing of cortical activities in multiple areas. The purpose of the present study was to investigate the temporal 
relationship among several cortical activities following ES using magnetoencephalography 

METHODS 

The experiments were performed in 13 healthy male volunteers, aged 28-42 (mean 34) years. To produce a noxious stimulus, ES was used. By 
pressing a pushpin-like needle electrode with a needle tip of 0.2 mm against the skin gently, the needle tip was inserted adjacent to the nerve endings 
of the thin myelinated fibers in the epidermis. A surface electrode 1 cm in diameter was placed on the skin at a distance of 4 cm from the needle 
electrode as the anode. The electric stimulus was a current constant square wave pulse delivered at a random interval of 0.1-0.3 Hz. The stimulus 
duration was 0.5 ms. The current intensity was 0.1-0.3 mA, which produced a definite pricking sensation. 

The evoked magnetic fields were measured with dual 37-channel axial-type first order magnetometers (Magnes, Biomagnetic Technologies, San 
Diego, CA). The two probes were centered on the C3 and C4 positions as based on the International 10/20 System. The magnetic fields were 
recorded with a filter of 0.1-200 Hz with a sampling rate of 2083 Hz. The analysis window was 100 ms before and 400 ms after the stimulus, and the 
prestimulus period was used as the direct current baseline. A hundred responses were collected and averaged. Sources of the recorded magnetic fields 
were determined using brain electric source analysis (BESA, NeuroScan, Mclean, VA). 

RESULTS 

A consistent magnetic field (termed 1M) was identified following ES in both hemispheres. Its peak latency was 149 ms for the hemisphere 
contralateral to the stimulated side and 166 ms for the ipsilateral hemisphere. Figure 1 shows recorded waveforms and results of the multi-dipole 
analysis in a representative case. The earliest activity was in the contralateral SI, which appeared slightly earlier than the onset of 1M. The ES-evoked 
SI activity consisted of three early components peaking at 88, 98 and 109 ms directed posteriorly, anteriorly and posteriorly, respectively, and a later 
component peaking at 151 ms directed posteriorly. The SI source was located in the anterior crown of the postcentral gyrus, probably area 1. At a 
latency period around 1M, SII and insula were activated. The onset latency of SII activity was 109 ms, which was later by 28 ms than that of the SI 
activity. The onset latency of the insular activity was slightly (7 ms) earlier than that of SII activity but the difference did not reach a significant level. 
The peak latency of SII and insular activities were almost the same as the peak of 1M in both hemispheres, therefore SII, insula and late SI activity 
contributed to shape 1M. To explain the magnetic fields later than 1M, sources in the anterior cingulate cortex (ACC) and anterior medial temporal 
area (MT) were necessary. Both activities started at around 150 ms and peaked 200 and 300 ms, that were clearly later than SI, SII and insula 
activities. 

DISCUSSION 

The present finding that ES-evoked SI activity clearly preceded SII activity suggested a serial processing in SI and SII, which is consistent with 
anatomical findings in monkeys that SI receives projections from the thalamus and, in turn, sends projections to SII. Although the existence of sparse 
projections from the thalamus to SII in monkeys is also known [Friedman, 1986], electrophysiological studies have shown that this thalamic inputs to 
SII in simian primates cannot activate SII neurons independently of SI [Pons, 1987], supporting the serial mode of processing. Probably, laser-evoked 
SI activities paralleled with SII activity in previous MEG studies [Ploner, 1999][Kanda, 2000] correspond to our late SI activity. 

The present result that SII activity appeared slightly later than insular activity indicated an origin other than SII for activation of the insula. Since 
the dorsal part of the anterior/mid insula is specifically activated by thermal and noxious stimuli in both fMRI [Schnitzler, 2000] and MEG [Inui, 
2003] studies, it seems possible that the dorsal insula receives inputs from modality-specific neurons in the thalamus. For example, Craig [Craig, 
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1994] demonstrated a very high concentration of pain- and thermo-specific 
neurons in the posterior part of the ventral medial thalamic nucleus (VMpo), 
which has dense lamina I spinothalamic terminations. VMpo projects to the 
insula and in fact, dorsal anterior part of the insula in monkeys contains a 
concentrated number of nociceptive neurons. The insula has been shown to 
project to linrbic structures including amygdaloid complex and cingulate 
cortex. Our results showed that onset latencies of the MT and ACC 
approximately corresponded to the peak latency of insular activity, so it is 
possible that both MT and ACC are driven by insula. Therefore, there are two 
pathways: thalamus (VPL)-SI-SII and thalamus (VMpo)-insula-MT and ACC. 
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Fig. 1. Temporal profile of cortical activities following epidermal electrical 
stimulation (ES) applied to the left hand. The upper three traces are 
superimposed waveforms recorded from 37 channels in both hemispheres and 
evoked potentials recorded in Cz. The lower seven traces are temporal profiles 
of each cortical activity. Filled circles indicate a group of early SI activities. 

SI, primary somatosensory cortex; SII, secondary somatosensory cortex; MT, 
anterior medial temporal area. 
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EEG and MEG responses following stimulation of unmyelinated C fibers 

Qiu, Y., Inui, K., Tran, T.D., Wang X., and Kakigi, R. 

Department of Integrative Physiology, National Institute for Physiological Sciences, Okazaki, Japan 


ABSTRACT 

Cerebral processing of first pain, associated with A-delta fibers, has been studied intensively, but the cerebral processing associated with C-fibers, 
relating to second pain, remains to be investigated. This is the first study to clarify the primary cortical processing of second pain by 
magnetoencephalography (MEG). We selective activated C-fibers by the stimulation of a tiny area of skin with a CO 2 laser. 

As for the primary component (1M), in the hemisphere contralateral to the stimulation, two regions in the hand area of the primary somatosensory 
cortex (SI) and secondary somatosensory cortex (Sll)-insular were activated. The onset and peak latency of the two sources in SI and Sll-insular were 
not significantly different. In the hemisphere ipsilateral to the stimulation, only one source was estimated in Sll-insular, and its peak latency was 
significantly longer than that of the Sll-insular source in the contralateral hemisphere, probably through corpus callosum. Our findings suggest that 
parallel activation of SI and Sll-insular contralateral to the stimulation represents the first step in the cortical processing of C-fiber-related activities. 
In addition to SI and Sll-insular, cingulate cortex and medial temporal area (MT) around amygdala and hippocampus in bilateral hemispheres were 
also activated for the subsequent component, 2M. 

All components of EEG and MEG responses were significantly reduced in amplitude during distraction and diminished during sleep, particularly 
2M component. These findings indicate that these regions are related to the cognitive aspect of second pain perception, particularly activities in 
cingulate cortex. 

KEY WORDS 

Magnetoencephalography (MEG); Pain; C fiber; Attention; Electroencephalogram (EEG) 

INTRODUCTION 

There are two systems for nociceptive perception, one ascending through A-delta fibers relating to first or sharp pain, and one ascending through 
C-fibers relating to second or burning pain. A C0 2 laser beam is frequently used to record pain-related brain potentials (laser evoked potentials, LEP) 
and magnetic fields (LEF) in humans, since it activates specifically nociceptors without stimulating mechanoreceptors relating to tactile sensation. 
Recently, [Bragard, 1996] [Opsommer, 1999] reported a noninvasive and simple method for recording ultra-late LEP by selective activation of C 
afferent sensory terminals in the skin using a C0 2 laser to stimulate a tiny surface area. The physiological background of this method is that the C 
afferent sensory terminals in the skin have a higher density and lower activation threshold than the A-delta terminals. We have recorded a clear ultra- 
late LEP by modifying the method and reported that the conduction velocity of the peripheral nerve and spinal cord following this specific 
stimulation was approximately 1-4 m/s, which is within the range for unmyelinated fibers [Tran, 2001] [Tran, 2002a] [Qiu, 2001] [Kakigi, 2003]. We 
reported that 3-dipole model of SI and bilateral SII sources could explain the 1M component of MEG responses following stimulation of C-fibers 
[Tran, 2002b], and that cingulate cortex and medial temporal area (MT) around amygdala and hippocampus in bilateral hemisphere were also 
activated for the subsequent component, 2M [Qiu, 2004], and also clarified the mechanisms underlying the effects of attention/distraction on second 
pain perception in detail by using EEG [Qiu, 2002] and MEG [Qiu, 2004], 

METHODS 

Thirteen healthy male volunteers participated in this study. For recording the ultra-late LEP and LEF by selective stimulation of C-fiber, we used 
a thin aluminum plate reported in detail in our previous papers [Tran, 2001] [Qiu, 2003], The stimulus intensity reaching the skin was 5-12 W with a 
20 ms duration, which caused pressure, touch, or slight burning pain (or second pain). The ultra-late LEF were recorded with two probes (BTi, San 
Diego, CA) placing at C3 and C4, with a 0.1-30 Hz bandpass filter, and digitized at a sampling rate of 1041.7 Hz. We used the brain electric source 
analysis (BESA) software package made by Scherg (BESA2000, NeuroScan Inc., McLean, VA) for the computation of theoretical source generators 
in a 3-layer spherical head model. The ultra-late LEP was also recorded with three exploring electrodes that were placed at Cz (vertex), C3 (around 
the hand area of SI in the left hemisphere), and C4 (around the hand area of SI in the right hemisphere) referred to linked earlobes (A1+A2). Ultra- 
late LEP were recorded under five different sets of conditions, (i) Control condition in which subjects did not have any task, (ii) Attention in which 
the subjects were instructed to mentally count the number of laser pulses paying close attention to the stimuli, (iii) Distraction in which the subjects 
were asked to calculate, (iv) Drowsiness (stage 1 in sleep) and (v) Stage 2 in sleep. Ultra-late LEF was recorded under two different conditions 
between Control and Distraction. 

RESULTS 

1. Effects of attention, distraction and sleep on ultra-late LEP: 

The ultra-late LEP, small negative component (Nl) and a large 
positive component (PI), were recorded. Since the peak amplitude of 
Nl-Pl was largest at the Cz in all subjects in the control condition, 
we mainly analyzed the results recorded there. Results of ANOVA 
indicated a significant relationship between the PI amplitude and 
attention level (P=0.0011). The Pi amplitude was slightly increased 
during the attention task (P=0.19) but significantly decreased during the distraction task (P=0.0068), relative to the control by post-hoc comparison. 
The PI amplitude was significantly larger during the attention task than the distraction task (P=0.0003) (left side in Fig. 1). The Pi amplitude was 
much decreased in 5 subjects, and was absent in the other 5 subjects during drowsiness. The mean amplitude decrease during drowsiness was 
significantly relative to the control condition (P=0.0044, n=5; P=0.0001, n=10) (right side in Fig. 1). 



Stimulation Stimulation 


O 0.5 1.0 1.5 sec 0 0.5 1.0 1.5 sec 


49 






S11 


significantly longer in latency in all the subjects, and the 
inter-hemispheric difference in mean latencies were 21±26 
ms (P<0.01) and 32±26 ms (P< 0.002) (Fig. 2). The RMS of 
1M and 2M was significantly smaller than in the Control, 

P<0.01 for 1M and P< 0.005 for 2M, in the contralateral 
hemisphere (Fig. 2). In the ipsilateral hemisphere, the RMS of 
1M and 2M was significantly smaller than the Control value, 

P<0.005 for 1M and P<0.002 for 2M (Fig. 2). 

Results using BESA: We made a 6-source model by 
locating Source 1 in SI of the contralateral hemisphere, 

Source 2 and 3 in SH-insula of bilateral hemispheres, Source 
4 in the cingulate cortex and Source 5 and 6 in MT of 
bilateral hemispheres (Fig. 2). The GOF value for this 6- 
dipole model in the Control was 92.9±1.1%. We tried to make 
other models by adding or changing sources, but this 6-dipole 
model was consistently most appropriate in all the subjects. 

Change of each source during Distraction: The current 
strength of SI and SII was significantly smaller in the 
Distraction than Control condition, P< 0.05 for SI and P<0.01 
for Sll-insula (Fig. 2). The current strength of the cingulate 
cortex and MT was significantly smaller in Control than in Distraction, P<0.001 for the cingulate and P< 0.05 for the MT (Fig. 2). The Sll-insula and 
cingulate were more effectively attenuated than SI and MT, particularly activities in cingulate cortex. 

DISCUSSION 

A painful laser stimulus activates both A-delta fibers and C-fibers, but the activation of A-delta fibers suppresses that of C-fibers. However, the 
cortical representation of the activation of unmyelinated C-fibers should play an important role in pain perception, relating to second pain and also 
visceral pain such as cancer pain. Our new technique, based on reports from a Belgian group [Bragard, 1996], enabled us to record clear ultra-late 
LEP and LEF. 

Since the subjects felt touch, pressure or a slight burning pain, we have to consider one fundamental and important question; Does this stimulus 
really relate to second pain? In conclusion, we believe the C-fiber responses in our study represent the second pain-related cerebral activation. There 
are 3 types of C-fibers, (1) polymodal or mechano-heat responsive C-fibers, (2) warm or heat responsive C-fibers, and (3) silent or mechano- and 
heat-insensitive C-fibers, and second pain is transmitted by (1). In normal subjects, C0 2 laser stimulation can activate (1) and (2). However, it is 
found that (2) have a higher heat threshold (48 °C) than (1) (40 °C). In the present study, we used a CO? laser to stimulate a tiny area of skin at weak 
intensity. Therefore, this stimulation probably activated polymodal C-fibers. 

The most important findings in this study is that we found activated regions following the stimulation of C-fibers and clarified their changes 
during the Distraction condition. Only one study has reported ultra-late LEF following C-fiber stimulation using the same stimulus methods as the 
present study (Tran, 2002b). They focused on the primary activities corresponding to our 1M component, and found that SI and SII were 
simultaneously activated in the contralateral hemisphere, and only SII was activated in the ipsilateral hemisphere. These findings were compatible 
with our results, and we added activities in the cingulate cortex, the MT, and probably the insula as mainly later activities in the present study. 

In conclusion, we confirmed that the cerebral response relating to signals ascending through C-fibers was much affected by the level of 
consciousness, possibly corresponding to second pain perception. Moreover, since the activities in the SI, Sll-insula, cingulate cortex and MT were 
significantly reduced when the subject was distracted from the stimuli, the level of attention should be strictly controlled when the ultra-late LEP and 
LEF recording technique is used for clinical purposes. 
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Neuroimaging of Pain 

K.D. Davis 

Toronto Western Research Institute & The University of Toronto, Canada 

Hemodynamic-based techniques (fMRI, PET) have revealed widespread areas of cortical and subcortical activations during application of painful 
stimuli. PET and fMRI indirectly relate to neuronal activity and cannot resolve the fine temporal details. But, the more direct high temporal 
resolution technique of MEG can be useful to study pain. The goal of this lecture is review the techniques used to examine pain and then consider 
analysis approaches and interpretations. Issues to be discussed with illustrative examples include (i) flavours of pain-related activations extracted 
using percept-related fMRI [1][2]; (ii) pain-related vs nociceptive-stimulus-related activations in hysterical anesthesia [3] and irritable bowel 
syndrome. 

(Supported by the CIHR, OMHF and KDD is a Canada Research Chair in Brain and Behaviour). 
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Psychophysiological state dependency of pain-relevant brain activity revealed by functional 

neuroimaging 

J. Lorenz 

Institute of Neurophysiology and Pathophysiology, University Clinic Hamburg-Eppendorf, Hamburg, Germany 

We investigated effects of expectation on intensity ratings and somatosensory evoked magnetic fields and electrical potentials following painful 
infrared laser stimuli in healthy subjects, using 31-channel magnetoencephalography (MEG) and 64-channel electroencephalography (EEG), 
respectively. Different auditory cues infomied subjects whether low or high laser intensities would follow. Twenty percent of trials were invalidly 
cued. In high correlation with perceived pain, the peak amplitude of MEG at a latency of 165 ms localized within SII and varied dependent on 
stimulus intensity and cue validity. The EEG maximum at a latency of 300 ms yielded a source in the posterior cingulate cortex (PCC) that also 
varied with stimulus intensity, but failed to show any cue validity effects. However, unlike SII activity PCC activity was generally greater during 
cued compared to control blocks without cues. In another experiment we used positron emission tomography (PET) to detect the brain anatomical 
basis of pain during inflammatory conditions. We applied a topical solution of capsaicin to the skin of healthy volunteers and used slowly-ramped 
contact heat stimuli. In one set of scans we applied a moderately painful stimulus upon the normal skin, in another set we applied the same subjective 
pain intensity with a normally warn stimulus on sensitised skin. Image subtractions comparing these two equally painful intensities reveal the 
specific activation of limbic forebrain structures such as anterior insula, perigenual anterior cingulate, ventral striatum, and prefrontal cortex during 
pain originating from sensitised skin. We conclude that top-down modulation of pain perception is governed by limbic brain structures and recruited 
according to the psychophysiological status of the body determined by anticipation, attentional state and pain type. 
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Advantages of MEG to Investigate Complex Epilepsy: An Overview 

N. Nakasato 

Department of Neurosurgery, Kohnan Hospital, Sendai, Japan 

Epileptic activity arises from a small part of the brain, extends locally, and propagates into remote areas within several tens of milliseconds, so 
localization of the primary source is clinically important. Scalp EEG and MEG can identify epileptic discharges only if the signal is distinguishable 
from the background activity. The diagnostic value of EEG and MEG depends on clear signal separation from the background brain activity and from 
remotely propagated or locally extended activity. Theoretically, MEG has higher spatial resolution than scalp EEG because radial current is not 
detected and inhomogeneous head conductivity causes little distortion. Several advantages of MEG have been shown clinically: MEG recognizes 
small tangential spikes overlooked by scalp EEG [1]; MEG can identify the primary epileptic region in patients with synchronized bilateral spikes 
[2]; MEG detects rhythmic spikes not detected by scalp EEG during seizures [3]; and benign rolandic spikes in childhood may be localized in the 
precentral cortex if surface negativity in gyral cortex, visible only in EEG, continuously extends to the MEG-visible fissural cortex facing the central, 
sylvian and interhemispheric fissures. 

[1] Park, H.M., Nakasato, N., et al. (in press) Comparison of magnetoencephalographic spikes with and without concurrent electroencephalographic 
spikes in extratemporal epilepsy. Tohoku J Exp Med. 

[2] Yu, H.Y., Nakasato, N., et al. (in press) Neuromagnetic separation of secondarily bilateral synchronized spike foci: report of three cases. J Clin 
Neurosci. 

[3] Iwasaki, M., Nakasato, N., et al. 2003. Focal magnetoencephalographic spikes in the superior temporal plane undetected by scalp EEG. J Clin 
Neurosci 10, 236-238. 


MSI yields an additional tool for successful presurgical evaluation of frontal lobe epilepsy. 

°P. Ossenblok, 'J.C. de Munck, 2 F.S.S. Leijten, 'A. de Jongh, °A. Colon, ° 3 P. Boon. 

“Epilepsy Centre Kempenhaeghe, Heeze, Netherlands; 'MEG Centre, VUmc, Amsterdam, Netherlands; 

^University Medical Centre, Utrecht, Netherlands; 3 University Hospital Gent, Gent, Belgium 

Diagnosis and surgical treatment of frontal lobe epilepsy (FLE) remains complex, because of poor electroclinical localization of this type of epilepsy. 
This study aimed at demonstrating that advanced source analysis of interictal MEG yields an additional tool for preoperative localization of FLE 
compared to EEG alone. 

The additional value of the interictal MEG versus EEG was assessed systematically for 21 patients with FLE, with regard to yield, sharpness and 
localizability of the interictal epileptifonn discharges (IEDs). Visual inspection showed that in general interictal MEG discharges were more 
abundantly present than in the EEG, while MEG spikes tended to be of shorter duration. Cluster analysis was used to classify these discharges on the 
basis of their spatial distribution followed by equivalent dipole analysis of the spike sub-averages, taking into account the realistic shape of the head. 
It was investigated systematically whether the grouping of spikes allowed for computation of spike sub-averages with enhanced signal-to-noise ratio 
and improved accuracy of source localisation. Furthermore, a model-based study predicted a higher spike yield in MEG versus EEG, because of the 
substantially higher SNR of the MEG in the frontal areas. However, insufficient SNR is not the only explanation for the failure of the localization of 
EEG spike clusters for almost twice as many patients as for the MEG. Independent source analysis of simultaneously recorded interictal MEG and 
EEG events yielded evidence that it is easier to differentiate the primary area of activity from the more widely distributed secondary propagation 
areas. Thus, in general MEG yielded a better “goodness of clinical fit“ compared, either to an epileptogenic lesion at MRI (n=18) and/or compared to 
intracranial recordings (n=10) than EEG. 
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Whole Head MEG and EEG in the Presurgical Evaluation of Epilepsy Patients: 

A Prospective Study 

S.Knake 1 ’ 2 , S. M. Stufflebeam 1 , Hideaki Shiraishi 1 , P. Ellen Grant 1 , Keiko Hara 1 , Thomas Witzel 1 , Valerie A. Carr 1 , Deirdre Foxe 1 , 
Evelina Busa 1 , Matti S. Hamalainen 1 , Seppo Ahlfors 1 , Susana Camposano 1 , Donald L. Schomer 3 , Edward B. Bromfield 4 , Barbara A. 
Dworetzky 4 , Joseph R. Madsen 5 , Reese G. Cosgrove 6 , Peter L. Black 4 , Blaise F. Bourgeois 5 , Pal G. Larsson 7 , Andrew J. Cole 6 , 

Elizabeth Thiele 6 , Anders M. Dale 1 and Eric Halgren 1 
'MGH/MIT/HMS Athinoula A. Martinos Center for Biomedical Imaging, Charlestown, MA 

Purpose: To evaluate the clinical impact of a combined 306-channel MEG / 70-EEG-recording on the presurgical evaluation of epilepsy patients. 
Methods: Sixtyseven patients were prospectively evaluated by simultaneously recorded MEG and EEG. MEG and EEG were evaluated by two 
independent reviewers. All patients were surgical candidates and had undergone Video-EEG-Monitoring. Two source modelling techniques were 
used: Equivalent current dipole fitting (ECD) and dynamic statistic parametric mapping (dSPM). 

Results: The overall sensitivity to detect interictal epileptiform discharges was 72% (48 / 67 pts) for MEG and 61 % for EEG (41 / 67 patients). The 
combined sensitivity was 75%. The simultaneous MEG/EEG investigation has induced a change in clinical management in 13 of 67 patients (19%). 
Conclusion: MEG and simultaneously acquired EEG are complementary and an enrichment in the presurgical evaluation. The MEG/EEG results 
impacted the clinical management in this selected group of the patients. The combined evaluation improved the selection of candidates for epilepsy 
surgery. Large prospective studies on the impact of MEG on the postsurgical outcome have to follow. 


MEG in Presurgical Evaluation 

Scheler G., Fischer M.J.M., Hummel C. Stefan, S. 

Epilepsy center and Neurocenter, Dept, of Neurology, University of Erlangen-Nuremberg 

Accurate assessment of the brain tissue generating epileptic activity is fundamental for epilepsy surgery. Several diagnostic methods participate in the 
process of presurgical evaluation, supplying infomiation on various morphological and functional aspects which is ultimately integrated into a final 
treatment decision. Magnetic source imaging (MSI), combining structural (magnetic resonance imaging, MRI) and functional 
(magnetoencephalography, MEG) data, has an increasing impact in presurgical epilepsy evaluation. MEG identified specific epileptic activity in 70% 
in a study of 455 epilepsy patients who underwent MSI investigations in our hospital. MSI localisation agreed with the surgically treated lobe in 89% 
of 131 patients, results for extratemporal cases being superior to those with temporal lobe surgery. A favourable outcome is correlated to a good 
match between resection volume and MEG results (presented at this meeting). Introducing a measure to quantify contribution of MSI to the general 
result of presurgical evaluation, the results showed that in more than one third, MSI had supplied additional information, and in 10% information 
crucial to decision making. Accuracy as well as contribution findings underline MSI ability to identify the irritative zone in temporal and 
extratemporal epilepsies. 
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Auditory Sensory Gating Deficit and Cortical Thickness in Schizophrenia 
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ABSTRACT 

Both an EEG P50 sensory gating deficit and abnormalities of the temporal lobe structure are considered characteristic of schizophrenia. The 
standard P50 sensory gating measure does not foster differential assessment of left- and right-hemisphere contributions, but its analogous MEG M50 
component may be used to measure gating of auditory source dipoles localizing to left- and right-hemisphere primary auditory cortex. The present 
study sought to detennine how sensory gating ratio may relate to cortical thickness at the site of the auditory dipole localization. A standard auditory 
paired-click paradigm was used during MEG data collection for patients (n=22) and normal controls (n=ll). Sensory gating ratios were determined 
by measuring the strength of the 50 ms response to the second click divided by the first click (S2/S1). Cortical thickness was assessed by two reliable 
raters using 3D sMRI. Results showed that: (1) the schizophrenia group had a P50 and left M50 sensory gating deficit relative to controls; (2) cortex 
in both hemispheres was thicker in the control group; (3) in schizophrenia, poorer left-hemisphere M50 sensory gating correlated with thinner left- 
hemisphere auditory cortical thickness; and (4) poorer right hemisphere M50 auditory sensory gating ratio correlated with thinner right-hemisphere 
auditory cortical thickness. The MEG-assessed hemisphere-specific auditory sensory gating ratio may be driven by this structural abnormality in 
auditory cortex. 

KEY WORDS 

P50 Sensory Gating, Cortical Thickness, Electroencephalography, Magnetoencephalography, EEG, MEG 

INTRODUCTION 

Schizophrenia is a disorder characterized by abnonnal information processing. Reduced gating of the P50 component of the event-related brain 
potential is reliably found in patients with schizophrenia [e.g., Adler et al., 1982] and is generally considered to be a neurophysiological marker of 
attentional impairment. Clinical studies of M50 gating have shown that gating predicts performance on neuropsychological tests of attention, working 
memory, general memory, and executive function in schizophrenia. Right-lateralized M50 sensory gating ratio is correlated with severity of negative 
symptoms in schizophrenia [Thoma et al., in press], whereas left-lateralized M50 gating ratio is related to positive symptoms [Irwin et al., 2003]. 

Temporal lobe structural abnormality in STG areas, thought to be generators of P50 and M50, and atypical asymmetry of planum temporale (a 
possible generator of P50) have been well documented in schizophrenia. There is abundant evidence that a deficit in the cortical structure on the 
superior bank of the left temporal lobe contributes to symptoms in schizophrenia and that neurophysiological abnormality is related to the reduction. 
Thus, it was hypothesized that, to the extent that patients with schizophrenia demonstrated impaired sensory gating, they would also have reduced 
cortical thickness measured at the localization of the source dipole on Heschl’s gyrus and/or planum temporale. That is, auditory cortical thickness, as 
an index of cortical pathogenesis, will be negatively correlated with sensory gating ratio. 

METHODS 

Subjects: Eleven normal control subjects and twenty-two patients with schizophrenia, matched on age and education, were solicited through ads 
in the local media. Following entry into the study, one normal control subject was excluded due to abnormalities found in his sMRI. Patients had a 
minimum of three months on their current medication and had not had an inpatient stay for at least one year, were current outpatients with chronic 
schizophrenia, had been hospitalized at least once, and were well characterized by their providers. Time since first hospitalization ranged from 10 to 
39 years (p = 23, sd = 9). All subjects were screened with the Structured Clinical Interview for DSM-IV Axis-I Disorders, Clinician Version (SCID- 
CV) to establish appropriate diagnosis, screen for co-morbid diagnoses, and screen for current substance abuse or history of drug dependence. 
Subjects were screened for history of neurological disorder or disease and serious medical illness. Also, control subjects were not included if they had 
a family history of a psychotic disorder in a first-degree relative. Screening procedures also included a "noise run," where subjects were placed in the 
MEG dewar to screen for metal-induced artifact. Prior to data collection, and after a complete description of the study to the subjects, written 
infonned consent was obtained. 

Neuroimaging Data Collection: MEG was recorded in a magnetically and electrically shielded room using a whole-brain, 122-channel 
biomagnetometer system [NeuroMag, Helsinki], At the start of each test session, subjects were fitted with an electrode cap to which three small coils 
were attached in order to provide specification of the position and orientation of the MEG sensors relative to the head. 150 uncontaminated tone-pair 
epochs were collected with the overall session lasting approximately 30 minutes. A bipolar oblique channel of electro-oculogram (EOG) was 
recorded simultaneously with MEG. Epochs were rejected if peak-to-peak signal strength exceeded 150 uv in any EOG channel and 3000 fT/cm in 
any MEG channel. MEG data were digitized at 500 Hz per channel for 1000 ms beginning 100 ms prior to SI. Separately, 3-D volumetric Tl- 
weighted sMRI was collected with a 1.5 T Picker Edge Imager. The pulse sequence was a Field Echo 3-D Sagittal (Picker) with the following 
parameters: TR =15 ms, TE = 4.4 ms, FOV = 256 mm, flip angle = 25 degrees, matrix 192 x 256, slice thickness =1.5 mm, no gap. 

Stimuli and Data Analysis: Stimuli were presented in pairs with 500 ms ISI and an intertrial interval that randomly varied. Stimuli were binaural 3 
ms clicks, created and delivered using NeuroStim software through Etymotic earphones placed in the subject's ear canal. Prior to data collection each 
subject's hearing threshold was detennined, and peak click intensity was set 30dB above threshold. In order to minimize motion-related artifact 
produced by the plastic tubes connected to the Etymotic earphones, the tubes were taped to the subject's face and ear. 

An equivalent current dipole model was used with the assumption that the neuronal sources were focal. Source locations were determined using 
NeuroMag (2000) software. To co-register MEG and MRI data, three anatomical landmarks (nasion and right and left preauriculars) were measured 
for each subject using the Probe Position Identification system [Polhemus] prior to data collection. The same three points were identified in the 
subject's sMRI, and a transfonnation matrix that involved rotation and translation between the MEG and sMRI coordinate systems was used. To 
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increase the reliability of the MEG-MRI co-registration, approximately 50 points on the scalp were digitized with the Polhemus system, in addition to 
the three landmarks, to ensure that all points were located on the MRI scalp surface. 

Using NeuroMag software, a 5 to 55 Hz bandpass filter and a -100 to -10 ms baseline adjustment were first applied to averaged MEG data. M50 
was defined as the first upward-oriented dipole occurring prior to Ml00, 40 to 80 ms poststimulus. Determination of the strength, location and peak 
latency of the M50 sources in left and right hemispheres was accomplished by fitting each dipole using subsets of 34 planar gradiometers over each 
temporal lobe. For modeling SI M50, 5 ms of data before and 5 ms of data following the M50 peak were selected. Equivalent current dipoles (ECDs) 
were then determined separately for each hemisphere using NeuroMag source localization routines. Only ECDs with goodness-of-fit values 
exceeding 75% for SI were accepted. Peak strength of the source over the 10 ms period was then determined. S2 M50 was identified using a 
procedure in which the location of the S2 dipole was assumed to be the same as that of the SI dipole, and the strength of S2 M50 was determined for 
the 10 ms period surrounding the peak latency obtained for SI. M50 sensory gating for each hemisphere was expressed as a ratio similar to that of 
P50: S2 dipole peak strength divided by SI dipole peak strength. 

Cortical Thicbiess Measure: In all cases, the 50 ms magnetic source dipole localization fell on posterior Heschl’s gyrus, on planum temporale, or 
within an expected 5 mm margin of error. NeuroMag MEG-MRI Integration Software [NeuroMag, 2001], designed to interact with dipole 
quantification and localization software, allowed linear measures of cortical thickness at the site of dipole localization. Two raters, blind to subject 
diagnosis and trained on NeuroMag source localization software, first achieved a 3-coordinate localization for the dipolar source (x, y, and z axes). 
The next task was to establish inter-rater reliability for cortical thickness measures. A tangent to the curve of the surface of the cortex at that 
localization was then approximated, and cortical thickness was measured as a linear distance perpendicular to that line. In cases where a fold in the 
cortex appeared to obscure the thickness, a similar procedure was used on a nearby area in which the thickness of cortex was more apparent. A level 
of .90 was set a priori as an acceptable level of inter-rater reliability. 

RESULTS 


Inter-rater reliability (intra-class correlation) for cortical thickness measures was.92 (p<.0001) and .91 (p<.0001) for left and right hemispheres, 
respectively. Consistent with our previous studies, patients had higher M50 sensory gating ratios in left (ANOVA F(l,31)=6.63, p=.02) but not right 
(F(l,31)=.006, p=.94) hemispheres. Patients had thinner cortex in both left (F(l,31)=90.68, p<001) and right (F(l,31)=49.58, p<.001) hemispheres at 
the site of dipole localization. To compare lateralized gating ratios with cortical thickness, left- and right-hemisphere M50 gating ratios were 
regressed on left STG cortical thickness and then on right STG cortical thickness. 

The prediction of left STG cortical thickness for the schizophrenia group accounted for 50% of the variance (F(2, 21)=9.32, p=.002). Evaluating 
each predictor for the unique variance it contributed, there was a negative correlation between left (p=.001; see Figure) but not right (p=.23) M50 
gating ratio and left STG cortical thickness. The prediction of right STG cortical thickness for the schizophrenia group accounted for 40% of the 
variance (F(2, 21)=6.24, p=.008). Evaluating each predictor for the unique variance it contributed, there was a negative correlation between right 
(p=.01; see Figure) but not left (p=.10) M50 gating ratio and right STG cortical thickness. 



Left Hemisphere M50 Auditory Gating Ratio 



Right Hemsiphere M50 Auditory Sensory Gating Ratio 


The prediction of left STG 
cortical thickness for the control group 
was non-significant. The prediction of 
right STG cortical thickness for the 
control group accounted for 71% of the 
variance (F(2, 9)=8.71, p=.01). 

Evaluating each predictor for the 
unique variance it contributed, there 
was a negative correlation between 
right (p=.01) but not left (p=.09) M50 
gating ratio and STG cortical thickness. 

Left and right STG thickness 
was correlated in patients (r = .71 ,P< 
.001) but only marginally so in controls 
(r = .51, p = .10), though there was no 
difference between these correlations 
(P=.67). 


DISCUSSION 

Patients with schizophrenia showed bilateral reductions in cortical thickness at auditory dipole localizations associated with M50 auditory 
sensory gating, a left-lateralized deficit in M50 sensory gating, and a hemisphere specific negative correlation between M50 gating ratio and cortical 
thickness. This correlation is reminiscent of McCarley et al. (2002) showing left-lateralized impairment in P300 in schizophrenia, bilateral reductions 
of gray matter volume on posterior STG, and a correlation between P300 and left cortical volume. No relationship between gating and thickness of 
cortex was expected in normal control subjects, and it was assumed that healthy cortex would be associated with unimpaired gating. This hypothesis 
held true for the left hemisphere, an unexpected negative correlation was found in the right hemisphere. This unlikely result was not predicted, and at 
the very least should be replicated prior to being interpreted. 

Thoma et al. (2003) linked impairment in auditory sensory gating with the function of a fronto-temporal network. Hanlon et al. (under review) 
showed that gating impairment is first seen at 50 ms in the left hemisphere and by 100 ms is seen bilaterally, suggesting that a deficit that begins very 
early in the left hemisphere rapidly spreads to become a more generalized deficit. The authors emphasized the role of a network in transmitting 
pathological information processing both temporally and spatially. It is tempting to conclude that the substrate for the overall gating impairment is in 
STG and that the signal measured there simply reflects abnormality of cortical thickness. Present analyses should be interpreted with caution, as the 
results are preliminary, and the correlations found between gross cortical structure (sMRI) and impaired cortical function (MEG) in schizophrenia 
may both reflect a more distal developmental cause. 
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Spectro-temporal neuromagnetic and hemodynamic correlates of impaired P50 suppression in 

schizophrenia 

K. Mathiak, H. Ackermann, A. Rapp, U. Klose, T.T.J. Kircher 
University of Tubingen, Germany 

Long standing observations on schizophrenia reveal several measures of stimulus adaptation to be impaired. The P50 suppression paradigm assesses 
the decrease of auditory evoked responses to a second stimulus S2 shortly (500 ms) after the first SI. Dysfunctions at the thalamus as well as cortical 
networks have been suggested as cause for impaired P50 suppression in schizophrenia. We studied 12 subjects with schizophrenia according to 
DSM-IV and a group of matched controls with P50 suppression paradigms in fMRI and whole-head MEG. Spectro-temporal analysis of evoked and 
induced neuromagnetic fields over the left and the right auditory cortex revealed reduced overall activity and a reduction of energy in the beta 
frequency range at later frequency ranges most pronounced over the right hemisphere. Moreover, the BOLD measure of suppression correlated most 
highly at the right hemisphere with its neuromagnetic equivalent across all subjects (cc = 0.5). Finally, coherences between the left and the right 
hemisphere of induced oscillations were significantly reduced in the patient group. The study reproduced stimulus suppression deficits in 
schizophrenia with localized BOLD responses and in neuromagnetic oscillation at the left and the right auditory cortex. Moreover, both measures 
were found to be linked by cross-subject regression. The data support models of impaired local network functions across primary and higher order 
auditory cortex in schizophrenia. These disturbances might contribute to the emergence of auditory hallucinations and formal thought disorders. 
Acknowledgments: Supported by DFG SFB 550/BI. 
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Inherited auditory-cortical abnormalities in twin pairs discordant for schizophrenia 
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ABSTRACT 

Information on the inheritance of neurophysiological deficits might help to elucidate the molecular genetic basis of schizophrenia. Here, we used 
MEG and EEG to determine inherited vs. state dependent auditory-cortical deficiencies in schizophrenia. Mono- and dizygotic (MZ and DZ) twin 
pairs discordant for schizophrenia, recruited from a total population cohort, were compared with a demographically balanced sample of healthy 
control twin pairs from Finland. EEG/MEG responses were measured to frequent standard and “rare” deviant tones. Neurophysiological measures 
were regressed against the degree of genetic relatedness to patients with schizophrenia. As expected, the EEG responses P50, Nl, and mismatch 
negativity (MMN), as well as the source amplitude of the MEG response P50m, were reduced in the schizophrenic patients. P50/P50m and N1 were 
significantly decreased also in their unaffected co-twins, as compared to the controls. Interestingly, in the unaffected subjects, the decrease of P50 and 
Nl amplitude correlated significantly with the degree of genetic resemblance to the schizophrenic subjects. Genetic modeling of variance 
components suggested that the P50 and N1 amplitudes have high heritability. In sum, schizophrenia may be associated with inherited deficits in early 
cortical auditory processing, reflected by the decreased P50 and Nl amplitudes. The MMN changes might reflect mainly disease-specific changes. 

KEY WORDS 

Attention; auditory; electroencephalography (EEG); magnetoencephalography (MEG), mismatch negativity (MMN), schizophrenia 

INTRODUCTION 


Despite the well-known genetic contribution [Cannon et al. 1998], the molecular genetic basis of schizophrenia still remains elusive. Auditory 
event-related potentials and magnetic fields (ERP and ERF, respectively) could help to identify endophenotypes, processes mediating between 
genotype and phenotype in schizophrenia. Alterations of the early component P50 generated at ~50 ms post stimulus, specifically the P50 gating 
deficits, are thought to reflect genetic predisposition to the disorder [Siegel et al. 1984], The Nl component and particularly the subsequent mismatch 
negativity (MMN), which peaks 100-200 ms after any detectable change in sound stimulation, are also reduced in schizophrenia [Javitt 2000]. Two 
familial studies suggested heritability of MMN deficits in schizophrenia [lessen et al. 2001, Michie et al. 2002], but contradictory findings also exist 
[Bramon et al. 2004, Salisbuiy et al. 2002], We studied whether abnonnalities in the auditory P50/P50m, Nl/Nlm, and MMN/MMNm covary with 
the degree of genetic relatedness to schizophrenics, in unaffected MZ and DZ co-twins of schizophrenic patients. 

METHODS 


MZ (10 pairs, 4 female, age 53 [47-62] years) and DZ (13 pairs, 7 female, age 51 [47-58] years) twin pairs discordant for schizophrenia (DSM- 
IV) and demographically balanced MZ (9 pairs, 4 female, age 51 [45-58] years) and DZ control pairs (11 pairs, 6 female, age 51 [47-58] years) were 
recruited from a cohort of all the same-sex twins born in Finland 1940-1957 [Kaprio et al. 2002], The institutional review boards of UCLA and the 
National Public Health Institute of Finland approved the protocols, and all subjects signed informed consent forms. The stimuli were 600-Hz standard 
pure tones (duration 50 ms, 5 ms rise/fall time), randomly replaced by 670-Hz deviants (P= 0.2), presented with a 450-ms inter-stimulus interval via a 
plastic tube to the right ear, 60dB over individual hearing threshold. In an additional block, backward-masking tones (randomly 300, 400, 900, or 
1000 Hz) followed each tone after a 125 ms gap. The subjects ignored the tones and concentrated on a silent video film. EEG and MEG were 
measured in a magnetically shielded room (Euroshield Ltd., Eura, Finland), the first 27 twin pairs with a 122-channel MEG device (passband 0.03- 
100 Hz, sampling rate 394 Hz) and the remaining 16 pairs with a 306-channel MEG device (passband 0.03-200 Hz, sampling rate 800 Hz) (Elekta 
Neuromag Ltd., Finland) [Ahonen et al. 1993]. Nose-referenced 64-channel (60 for the last 16 pairs) EEQ with bipolar vertical and horizontal EOG 
(passband 0.5-30 Hz), was recorded with a MEG-compatible electrode cap and amplifier [Virtanen et al. 1996, Virtanen et al. 1997,]. The locations 
of marker coils in relation to the subjects’ head were determined before measurements (Isotrak 3D-digitizer, Polhemus, Colchester, VT) to estimate 
the subjects’ head position in relation to the MEG instrument before each stimulus block. Off-line bandpass was 5—40 Hz for P50/P50m, 1-30 Hz for 
Nl/Nlm and MMN, and 2-15 Hz for MMNm. At least 600 artifact-free standard and 150 deviant epochs (750 ms, 150-ms baseline) were averaged. 
The P50 and Nl peak amplitude was determined from the average of electrodes FC1, FCz, FC2, Cl, Cz, and C2. The 100-200 ms post-stimulus 
average MMN (deviant minus standard ERP) amplitude was detennined from 6 frontocentral electrodes. The ERF peak latencies were determined 
from the peak total-field power, and the amplitudes by least-squares ECD fitting (spherical head model, origin {x, y, z}={0, 0, 45 nun}), based on a 
left temporal subset of gradiometer channels. The group differences in ERP/ERF were analyzed by using a mixed-model ANCOVA (age, gender, and 
recording system as covariates, twin pair as a random variable) with contrasts. The ERP/ERF were regressed against the degree of genetic relatedness 
to schizophrenic patients, across the unaffected subjects, using a mixed model with the twin pair as a random variable. The heritability of ERP was 
analyzed with genetic modeling of variance components, with age, gender, and the disease status as covariates (Mx [Neale 1997]). 

RESULTS 


According to the mixed-model contrasts, the P50 amplitude was significantly larger in the controls than schizophrenic patients (t (18.4)=3.6, 
p<0.01) or their co-twins (t (18.4)=3.7, /><0.01) (Fig. 1). Nl was significantly larger in the controls than schizophrenic patients (t (27.3)=3.1, 
/><0.01) and their co-twins (t (27.3)=3.2, j?<0.01). There were no P50 and Nl differences between the schizophrenics and their co-twins. The MMN 
amplitude was smaller in the schizophrenic patients than in their co-twins (t (31.1)=3.5,/><0.01) and controls (t (21.2)=3.2,/><0.01), but the co-twins 
and controls had no differences (Fig. 2). The P50m source amplitude was significantly larger in the controls than schizophrenics (t (17.7)=3.0, 
/><0.01), and there was a similar trend in the co-twins (p<0.10). P50m originated more posterior in the schizophrenic patients than controls (t (9.1)= 
-2.3, p<0.05), and there was a similar trend in the co-twins (p<0.10). The P50 (/?= -0.37,/><0.01) and Nl (/?=0.37, p<0.01) amplitude decrease 
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correlated with genetic similarity to schizophrenic patients (see also z-scores in Fig. 3). The genetic relatedness to schizophrenic patients tended to 
correlate, non-significantly, with the P50m source anterior-posterior location (J3= -0.35, /><0.10) and amplitude (Jh— 0.26, /><0.15). Genetic 
modeling of variance components showed high heritability of the P50 amplitude (h 2 = 0.55) and latency (/j 2 =0.59), and N1 amplitude (/j 2 =0.70). 



MZ discordant pair. 55 years 


DZ discordant pair. 57 years 


Figure 1. Single-subject scalp potential distributions at 
the P50 peak latency projected over a spherical head 
model in representative age-matched MZ/DZ twin pairs. 
P50 is reduced in the schizophrenics and their co-twins 



Figure 2. Grand-average deviant minus standard ERPs in the 
schizophrenic patients, their co-twins, and the controls. 
MMN was abolished in the schizophrenic patients, while no 
differences emerged between their co-twins and the controls. 



Controls DZ co-twins MZ co-twins 
Figure 3. Mean (SEM) z-scores of the 
ERP amplitudes in the healthy controls, 
and in the DZ and MZ co-twins of the 
schizophrenic patients. 


DISCUSSION 


Small amplitudes of the auditory P50 and Nl, shown to be highly heritable measures of brain function by our genetic modeling of variance 
components, correlated significantly with the unaffected subjects’ genetic relatedness to patients with schizophrenia. Further, there was a non¬ 
significant tendency toward analogous effects in the source location and amplitude of the MEG response P50m. Group-wise, the amplitude of P50 
and N1 was significantly decreased in both schizophrenic patients and their co-twins, and analogous non-significant trend was observed in the MEG 
responses P50m and Nlm. Previous studies have associated genetic predisposition to schizophrenia with P50 gating deficits [Siegel et al. 1984], A 
recent MEG study, however, suggests that the P50 gating deficits are explained by amplitude reduction in the response to the first stimulus in the pair 
[Blumenfeld et al. 2001]. Our results support the notion that the P50/P50m generation per se might be deficient in schizophrenics, further suggesting 
that this might be an inherited phenomenon. Supporting previous findings [Javitt 2000], the MMN amplitude was reduced in schizophrenic patients, 
but no evidence of genetic basis for this effect was found. This finding contrasts with some results [lessen et al. 2001, Michie et al. 2002], but can be 
supported by a very recent familial study [Bramon et al. 2004] and comparisons between first-episode and chronic schizophrenia patients [Salisbury 
et al. 2002], In conclusion, the decreased auditory P50 and Nl amplitude, and the changes in the P50m sources, might be related to genetic 
predisposition to schizophrenia. MMN in turn appears to reflect state-dependent changes in the schizophrenic brain. 
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A specific reduction in the sustained 40Hz magnetic field power in schizophrenia 

S. M. Stiifflebeam, T. Witzel, F. H. Lin, D.M. Foxe, K. Spencer, and R. C. McCarley 
HST/MIT/MGH Athinoula A. Martinos Center for Biomedical Imaging, Boston, MA, USA 

Purpose: To quantify and localize the neuromagnetic 40Hz sustained response in schizophrenia. 

Methods: Subjects : 5 medicated chronic schizophrenic subjects and 5 neurologically normal subjects were studied. Stimuli: Blocks of 500ms trains 
of 1msec auditory clicks were presented at 75dB HL, bilaterally. MEG: 306-channel VectorView; MRI: 3T fMRI using a sparse auditory imaging 
technique (sounds presented in 8 sec of silence). Analysis: Welch-method PSD using a 2-dipole ECD; Inflated cortex (FreeSurfer) constrained, 
wavelet-transformed minimum norm estimate with fMRI-weighting (Dale et al, Neuron 2000) to calculate current distributions on the cortical 
surface. 

Results: The 40Elz response was significantly reduced in the SZ subjects (P<0.05), but the 20 Elz and 30Elz response were not significantly different 
from the control subjects (P >0.7). fMRI demonstrates activation of the entire EleschTs gyrus with variable response in the planum temporally. There 
was a trend towards a correlation between the cortical thinning in the left posterior STG and the power reduction in the left 40EIz power, but did not 
reach statistical significance with our small pilot sample size. 

Conclusion: A specific reduction of the 40Hz sustained neuromagnetic response power in schizophrenic patients appears to originate in Heschl’s 
gyrus, most prominently on the left side. This supports findings in EEG [1] and suggest abnormal gamma processing the auditory cortex of 
schizophrenics. 

Supported by the MIND Institute. 

[1] Kwon et al. Gamma frequency-range abnormalities to auditory stimulation in schizophrenia. Arch Gen Psychiatry. 1999 Nov;56(ll): 1001-5. 
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Multimodality Studies of Simple Movement 

Mark Hallett, M.D. 

Human Motor Control Section, N1NDS, NIH, USA 

The combination of multiple physiological methods in the analysis of brain function reveals more than the individual measures, and can lead to better 
understanding. I will present two examples of the combination of transcranial magnetic stimulation (TMS), electroencephalography (EEG) and 
functional magnetic resonance imaging (fMRI). 

What happens in the contralateral motor cortex in the no-go task in a go/no-go task? EEG reveals activation similar to the go task. fMRI does not 
show any activity at all. TMS shows a reduction of excitability. Using TMS for testing of intracortical circuitry shows an increase in the short 
intracortical inhibition (SICI) without a change in short intracortical facilitation (SICF). These measures suggest that there is active inhibition in the 
motor cortex in the no-go task. It appears that EEG, in this circumstance, does not differentiate excitation from inhibition, and that fMRI does not 
detect inhibition effectively. 

What happens in the ipsilateral motor cortex during a simple unilateral movement? EEG (including power analysis, event-related desynchronization, 
ERD) shows some activity. fMRI shows only minimal activity. The motor evoked potential (MEP) produced by TMS in the contralateral hand is 
inhibited. Again, the interpretation is that involvement of the ipsilateral cortex is largely due to inhibition. 


Neural connectivity and its modulation by transcranial magnetic stimulation of the human frontal cortex. 

Tomas Paus 

Montreal Neurological Institute, McGill University, Montreal, Quebec, Canada. 

In the last two decades, two powerful tools for investigating brain mechanisms of cognition have emerged: functional neuroimaging and transcranial 
magnetic stimulation (TMS), the fonner capable of measuring and the latter of changing activity in the human brain. In this talk, I shall provide a 
brief overview of current trends in combining neuroimaging and TMS and the potential of this approach in studying brain plasticity. I will describe 
two sets of studies of TMS-induced changes in excitability and connectivity in the motor system, one based on the combination of TMS with positron 
emission tomography, and another on the combination of TMS with electroencephalography. 
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Combined EEG-TMS technique to study brain-behavior relations 

Alvaro Pascual-Leone, Gregor Thut 

Center for Non-lnvasive Brain Stimulation, Beth Israel Deaconess Medical Center and Harvard Medical School, Boston, MA 

Combined EEG-TMS techniques provide neuroscientists with a unique method to test hypothesis on functional connectivity, as well as on 
mechanisms of functional orchestration, reorganization, and plasticity. In addition, EEG monitoring during repetitive TMS, can be valuable to assess 
the safety of TMS and might be useful to monitor the modulation of cortical excitability by rTMS and adjust the stimulation parameters so as to 
achieve a given physiologic effects. However, EEG recording during TMS is technically challenging particularly because of the TMS-induced high- 
voltage artifacts. The use of slew-rate limited preamplifiers helps prevent saturation of the EEG system due to TMS and an artifact subtraction 
protocol allows isolation of the electrophysiological signals from residual TMS-induced contaminations. The reliability of the approach in the single¬ 
pulse TMS design, as well as during repetitive TMS can be illustrated with examples of two data sets. The effects of occipital single-pulse TMS on 
high-amplitude potentials evoked by visual checkerboard stimulation allows the characterization and reveals the duration of the impact of TMS on 
local cortical activity. The modulatory effects of rTMS on a neural network an be examined by measuring the impact of stimulation on low-amplitude 
steady-state oscillations induced by auditory click-trains. The main field of application of the present protocol for combined EEG-TMS studies is in 
cognitive neuroscience, complementing behavioral studies that use TMS to induce transient, ‘virtual lesions’. 


Linking changes in cortical excitability with changes in neuronal activity using a combined TMS- 

neuroimaging approach 

H.R. Siebner 

Department of Neurology, Christian-Albrechts-University, Kiel, Germany 
Neuroimage-Nord, Hamburg-Kiel-Liibeck, Germany 

Transcranial magnetic stimulation (TMS) can be used to quantify the excitability of the primary motor cortex (Ml) and, to some extent, of the visual 
cortex. For the Ml, cortical excitability is usually assessed further “downstream” in the motor system by measuring the transcranially evoked motor 
responses. In recent years, combined neuroimaging-rTMS studies have significantly extended the scope of currently available TMS techniques that 
probe cortical excitability because the combined neuroimaging-TMS approach can be applied to other cortical areas. Neuroimaging during TMS also 
provides a more direct probe of cortical excitability than conventional TMS techniques because BOLD signal, regional cerebral blood flow or 
regional cerebral metabolic rate of glucose are tightly coupled to regional neuronal activity. This talk will review a series of recent experiments in 
which PET or fMRI were combined with TMS of the Ml or dorsal premotor cortex. The results demonstrate that the combined TMS-neuroimaging 
approach provides a powerful means of linking changes in regional cortical excitability with changes in regional neuronal activity in the human brain. 
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